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PRKFACE 


Failure  analysis  of  fiber  reinforced  composite  itHterlals  Is  not  a  simple 
problem.  There  Is  a  wide  dimensional  ran;»e  In  which  analytical  failure 
theories  can  be  developed,  depending  on  the  viewing  scale.  Composite  loatar- 
lals  can  be  perceived  simply  as  a  class  of  homogeneous  elastic  media,  or  as  a 
highly  complex  structural  system.  There  are  a  variety  of  other  models 
proposed  In-between.  It  seems  Impossible  to  develop  an  adepuate  failure 
analysis  methodology  that  can  be  applied  over  the  entire  dimensional  spectrum 
of  composites  as  an  engineering  material. 

for  a  class  of  epoxy-based  structural  laminates,  failure  usually  begins 
with  mlccocracklng  in  the  resli  or  in  the  resln-flber  interfaces.  The  dimen¬ 
sions  of  the  mlcrocracks  are  on  the  order  of  the  reinforcing  fiber  llameter. 
During  the  course  of  loading,  some  of  the  mlcrocracks  may  coalesce  and  form 
what  Is  known  as  sub-lamlnate  matrix  cracks,  whose  can  he  orlers  of 

magnitude  larger  than  the  mlcrocracks.  Despite  such  large  sizes,  sub-l.amlna te 
matrix  cracks  generally  remain  localized  or  propagate  stably  until  reaching  a 
certain  late-staga  loading,  at  which  tine  n  rapid  g.rowth  of  the  'aatrl<  cracks 
occurs.  This  event  usually  renders  the  laminate  structurally  unsafe. 

The  ourpiase  of  thla  work  Is  to  liv*;stlgate  the  formation  -and  prnpavitlon 
mechanisms  of  some  Important  types  of  sub-l.aml  na te  matrix  cracks.  A  phennmon- 
ologlcal  failure  analysis  methodnlogv  Is  developel  within  t.’it>  Ilnenslnnnl 
bounds  la  which  these  matrix  c>'acks  a*"*  observed. 

It  should  be  mentlonel  that  composite  materials  possess  ninv  i\Iiue 
physical  characteristics  not  prevloijslv  enrount‘r?d  by  the  analyst^.  T'in;,  It 
Is  o'‘tan  necessary  tn  modify  the  existing  analysis  methods  for  each  uej’/ 
encountered  situations.  This  practice  miv  In  tira  cause  the  application  of 
the  analysis  methods  with  In-suf  f  iclent  scientific  uracl  slon  aal  m  i  themt  1 1 ‘■•a  1 


K?r, 
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rigor.  While  this  work  Is  not  an  exception  to  this  predicament,  It  Is  hoped 
that  It  will  stimulate  new  theoretical  research  In  the  field. 

The  author  Is  deeply  Indebted  to  his  colleaques,  former  and  present 
students;  through  whose  collaboration  much  of  the  results  presented  here  were 
obtained.  He  Is  especially  grateful  for  the  financial  supports  received  from 
the  H.  S.  Mr  Force  Office  of  Scientific  Research,  Wright  Aeronautical  Labora¬ 
tory  and  the  Naval  41r  Oevelopraent  Center  during  the  various  phases  of  his 
research. 

This  manuscript  was  prepared  while  the  author  was  on  a  sabbatical  leave 
from  Drexel  University  to  further  pursue  these  research  objectives  at  the 
Naval  Mr  Development  Center  under  a  program  sponsored  by  the  National 
Research  Council  from  'October  l,  1983  to  September  30,  1984.  Many  thanks  are 
due  to  Mr.  Lee  W.  Cause  of  NADC  for  numerous  helpful  technical  dlscussl'cns  on 
the  subject  of  damages  la  composites. 

Finally,  sincere  gratitude  Is  expressed  to  Mrs.  Shirley  Cheeseman  for  her 
patience  and  expert  typing  la  preparing  the  manuscript. 
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CHAPTER  I.  MATRIX  CRACKS  IN  COMPOSITE  LAMINATES 

One  rnalor  class  of  structural  composites  Is  fabricated  usln^  unMtrec- 
tlonal  plies  of  res  In-Impregnated  fibers  (prepreg  tape'l.  These  prepreg  piles 
may  be  stacked  In  a  certain  prescribed  orientation  sequence  to  achieve  the 
required  structural  properties  when  pressed  to  a  laminate.  Ideally,  the 
fibers  would  be  orientated  primarily  In  the  principal  loading  directions  with 
the  resin  taaterlal  serving  only  to  bind  the  fibers  and  the  various  olles 
together.  Experiences  have  shown,  however,  that  aaterlal  ‘'allure  la  this 
class  of  lamlnatas  seldom  begins  with  breaking  of  the  load-carrying  ‘'Ibers. 
Rather,  small  micro-cracks  first  occur  in  the  resin  or  In  the  resln-flher 
interfaces.  ’Inder  snstilied  loading,  these  mlcrocracks  can  form  In  increasing 
density  and  grow  larger  la  sice.  At  stress  critical  locations,  coalescence  of 
the  mlcrocracks  can  t-ake  place,  forming  what  is  known  as  matrix  cracks. 

Recause  of  the  lamlnata’s  Internal  reinforcement  structire,  matrix  cracks 
are  generally  Inhibited  from  monotonl:  propagation ,  but  they  ten.i  to  form  a 
cracking  pattern  w’llch  Is  characteristic  to  the  laminate  In  question. 

Experiments  on  laboratory  couoons  (see,  e.g.  ill)  have  established  that 
matrix  cracking  natterns  are  getierlc  to  the  ply  fiber  or  lent  atlons  and  stack¬ 
ing  sequence  of  the  laminate.  The  patterns  are  essential Iv  reiarodnclhle  la 
laminate  replicates  suhlected  to  similar  loading  con'll  tlons.  It  has  oft^n 
been  snggeste'i  that  the  basic  mechanisms  of  matrix  crackin'’  aa  i  the  assoclate'i 
characteristic  cracking  pattern  are  aneiaable  t>  some  I'orm  of  aiilvtlcal 
description  based  on  the  principles  of  mcctiaoi -rs  I  . 

From  the  viewpoint  of  'aaterlal  stran-gth,  a  crack  forned  In  flhro'is  'oiipii- 
sltes  signifies  that  the  limit  of  local  material  strength  has  Rees  ax.'''io.i. 
An  analysis  of  strength  shonll  ther;fore  he  cmlict^'i  at  thca  ni'roscooln  seal’ 
whera  the  composite  Is  dlstinctl/  a  mill  ti-oha-s  ’  iiterlal.  If  this  .'law  Is 
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I 

I 

I 

I 

i 

followed,  not  only  should  the  Individual  niechanlcal  properties  of  the  fleers 
the  matrix  and  the  fiber-matrix  Interfaces  be  Included  amon^  the  Innortsnt 
factors  determlnlnj’  the  streui^th  behavior,  but  local  inomslles  o'"  either 
geometrical  or  'natarlal  origins,  or  any  other  such  quantities  whose  Hmenslons 
are  comparable  to  the  fiber  diameter  should  also  be  considered.  One  malar 
difficulty  with  this  viewpoint  Is  that  material  anomalies  are  rendon  In 
nature,  rendering  fruitless  efforts  to  Include  their  Identity  and  e‘"fects  la 
any  deterministic  form  of  analysis. 

Yet,  mlcroflaws  do  exist  In  composites  and  they  remain  the  primary  cause 
of  matrix  cracking.  Their  presence,  or  at  least  their  effects,  must  In  some 
way  be  Included  In  the  analysts. 

4n  alternative  approach  Is  t'J  look  upon  the  matrix  cracking  process  as  a 
phenomenological  event.  Typically,  the  Individual  material  piles  are  regaried 
as  homogeneous  anisotropic  bodies.  In  which  the  actual  f Iber-ma t -I x  alcro- 

1 

I  structure  as  well  as  the  micro-defects  all  become  ladlstlngulshable.  Material 

I 

property  distinction  Is  ma fata  lied  between  plies  of  differing  fiber  orienta¬ 
tions.  Thus,  a  transition  Is  made  froii  a  nlcroscoptc  to  s  mcroscoiolc 
description,  with  the  requirement  that  Important  Factors  In  the  mlcroscinlc 
regime  be  retailed  In  the  macruscoplc  regime.  In  this  context,  t'le  orisenr*^ 
and  the  effects  of  material  mlcroflaws  on  matrix  cracking  nrs  represented  hy  a 
distribution  of  "effective"  material  flaws,  with  the  eFFectlve  Flaw  dlstrlhu- 
tlon  characterised  as  an  Inherent  property  of  the  miterlal  plv.  In  o’-l.jr  ti 
develop  a  methodology  based  oa  this  anproich,  hj\/eoer,  a  plivslcal 
Ini'  of  the  matrix  cracking  mechanisms  at  both  the  ifrroscoplc  and  f  ie  ii'r  >- 


scoolc  scale  Is  essential. 
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1 . I  Physical  Charicterlstlcs  of  Matrix  Cracks 

Because  of  the  unique  Internal  reinforcement  structure  of  the  consMerel 
class  of  laminates,  there  are  only  two  basic  forms  of  matrlK  cracklnq  when 
viewed  at  the  'sacroscale.  Suppose  that  the  stresses  la  each  ply  of  the 
laminate  ace  three  dimensional  [3];  then  there  are  three  In-ply  stresses  and 
three  Interply  stresses,  as  shown  In  Figure  l.l.  Under  certain  conditions, 
the  In-ply  stresses  may  Initiate  and  propagate  m-itrl'c  cracks  parallel  to  the 
fibers  In  that  ply,  known  commonly  as  tran3./erse  cracks,  Figure  1.?.  On  the 
other  hand,  the  Interply  stresses  may  Initiate  and  propa,qate  matrl<  cracks 
parallel  to  the  plane  of  the  ply,  or  In  the  Interface  between  two  piles.  The 
latter  l3  labeled  as  delamlna tlon.  Figure  1.3. 

The  moilnl  behavior  of  laminate  matrl<  cracklnq  clearly  depends  on  the 
relative  dominance  between  the  In-ply  stresses  and  the  Interply  stresses.  In 
many  instances,  the  state  of  the  ply  stresses  chany'es  as  local  cracks  form. 
Then,  a  combination  of  the  two  different  cracklnr'  forms  can  he  Induceil. 

'iul t Ipl Id ty  la  another  nnlque  character  of  ivntrl<  cracklny.  This  pheno¬ 
menon  la  aqaln  dne  to  the  laminates  Internal  'elnforcement.  Becanse  of  the 
filamentary  stnictire  In  each  ply  anil  tie  layerlny  structire  In  the  laminate, 
a  natrl>c  crack  cannot,  generally,  undi»r./o  nonotonlc  pronnyntlon.  The  crack  fa 
likely  to  be  arrested  or  blunted  at  aoiie  r  el  nf  nrceiient  Ixionliry.  This  crack- 
arrestln.q  mechanlam  permits  th''  applied  load  ti  Incrreiie  without  •.ansli<’  the 
crack  tn  propiynte.  The  resiit  Is  that  another  ivitrl<  crack  Is  lit  tint’d 
elsewhere  la  the  la'nlnate.  In  this  iinner,  a  1  in  l-dmcndcn t  niltlpT’  ^rick 
formation  process  la  ievelopel. 

Tnns,  It  is  cl>enr  that  the  ph/slcal  chnrioter  if  nitrl<  crio-  In'’  la 
Inflnenoel  atron.’!/  hy  the  laminate's  Internal  and  e<t’riil  •’I’omi' t  r  [ 's ,  Incln- 
lln;’  ply  fiber  orientations,  o1  /  thi  I'kiaess ’s ,  nl  /  stackin’  saii'-'ic’  an  I 
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laminate  shape  feature,  etc.  Consequently,  any  inayltlcal  approach  to  this 
problen  must  Include  Implicitly  or  explicitly  these  (geometrical  factors. 

1.2  Statistical  Homogeneity  and  Ply  Elasticity 

The  mathematical  transition  from  a  microscopic  tlescrlptlon  to  s  macro¬ 
scopic  tlescrlptlon  of  the  composite  laminate  Is  made  through  the  concent  of  i 
statistically  homof^eneous  body  [A].  This  concept  applies  to  flhrinns  compo¬ 
sites  In  general  and  to  the  nnl  ll  cectlonal  material  svste  ns  In  nartlciPar.  It 
Is  assumed  that  the  constltutlv/e  properties  of  a  unidirectional  olv  haying 
arbitrary  finite  volume  V  can  be  determined  from  the  same  ntcrlsl  having  a 
smaller  representative  volume  AV.  However,  In  order  for  the  determined 
properties  to  be  statistically  representative  of  those  In  volume  V,  the 
minimum  slae  of  aV  must  be  large  enough  to  Include  all  possible  variations  of 
the  composite  microstructure.  Including  fiber  packing  patterns,  presence  of 
microflaws,  etc.;  see  Figure  1.4. 

Within  this  context,  two  mutually  equivalent  approaches  can  be  taken  to 
characterize  the  constitutive  properties  of  the  unidirectional  ply.  The  first 
approach  Is  to  apply  the  methods  of  mlcromechanlcs  which  describe  the  mechani¬ 
cal  behavior  of  the  composite  of  volume  AV  In  terms  of  the  properties  of  the 
fiber  and  the  'aatrlx  phases  as  well  as  the  geometrical  configurations  of  the 
fiber-matrix  microstructure.  An  average  is  then  taken  for  the  stresses  and 
the  strains  In  AV,  whereby  a  set  of  effective  constitutive  properties  Is 
obtained.  By  virtue  of  the  statistically  homogeneous  bodv  pos tula t Ion ,  these 
effective  properties  are  regarded  as  the  same  for  the  entire  material  ?yst-’n 
of  any  arbitrarily  large  volume. 

The  other  approach  Is  to  deter  nine  phenomenologically  the  effective 
properties  of  the  material  ply  In  a  large  volume  hv  llrect  experimental 
Then,  ov.;lag  to  the  stati  cally  hiamgeiaeoiis  hodv  postiVite,  the s » 
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effective  properties  are  assumed  the  same  for  any  arbitrarily  small  volume  of 
the  material,  Including  those  In  aV. 

Traditionally,  the  effective  elastic  properties  determined  by  the  two 
different  approaches  have  been  treated  as  mutually  equivalent.  But,  extension 
of  this  equivalence  to  material  strength  oropertles,  and  especially  Its  impli¬ 
cations  on  the  mathematical  failure  criteria  have  not  been  critically  Investi¬ 
gated. 

Consider  for  example  a  typical  graphite-epoxy  system.  ’^‘Igure  1.5  shows  a 
cross-sectional  view  of  the  fiber-matrix  microstructure  enlarged  500x.  The 
thickness  of  the  Individual  material  ply  Is  In  the  order  of  125  ym,  whereas 
the  fiber  diameter  is  about  B  ym.  From  the  statistically  homogeneous  body 
concept  discussed  above,  the  minimum  dimension  of  a  representative  volume  aV 
should  contain  several  fiber  cross-sections  (say,  10').  If  so,  the  linear 
dimension  of  AV  Is  about  20Z  the  thickness  of  the  material  ply.  It  Is  thus 
seen  that  the  physical  size  of  AV  Is  quite  large  compared  to  the  ply  thick¬ 
ness.  This  raises  the  serious  question  as  to  whether  the  ply  homogeneity 
assumption  can  adequately  reflect  the  discrete  field  (microscopic)  effects 
when  a  crack-llke  failure  propagates  Inside  the  thickness  of  the  ply,  such  as 
In  the  transverse  cracks  shown  In  Figure  l.l.  Some  of  the  questions  regarding 
the  ply  homogeneity  assumption  and  the  use  of  the  fracture  mechanics  approach 
do,  Indeed,  pose  conceptual  difficulties. 

Consider  the  theory  of  ply  elasticity  [5],  an  Important  derivative  of  the 
statistical  homogeneous  body  concept.  The  theory  assumes  linear  elasticity 
for  the  !iaterlal  ply,  In  addition  to  the  homogeneity  assumption.  Specific¬ 
ally,  the  generalized  Hookes  Law  Is  assumed  to  govern  the  ply  constitutive 
relations: 

01  Cij  ej;l,1=l.6  (l.l) 
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wher'3  the  stiffness  constants  Cjj  are  determined  either  from  a  mlcromechanlcs 
analysis  method  for  a  representative  volume  AV  (see,  e.g.  [4]),  or  by  a  set  of 
experimental  test  methods  for  a  large  volume  V  of  the  ply  material  (see,  e.g. 
(6]).  In  either  manner,  the  determined  constants  C^j  are  assumed  a  uniform 
property  for  the  unidirectional  material  system. 

For  laminates  fabricated  by  stacking  piles  of  different  fiber  orienta¬ 
tions,  the  analysis  allows  the  ply  material  constants  to  change  abruptly 
across  the  ply  Interfaces.  Then,  assuming  a  perfect  bonding  at  the  Inter¬ 
faces,  an  exact  analysis  of  the  laminate  stress  field  based  on  the  ply  elasti¬ 
city  theory  can  be  performed  [31.  Such  an  analysis  provides  the  stress 
distribution  In  each  of  the  piles  from  which  a  suitable  ply  failure  criterion 
Is  applied  to  determine  whether  or  not  any  of  the  piles  Is  failed  by  the 
applied  load.  Typically,  the  ply  failure  criterion  Is  derived  from  the 
conventional  point-stress  concept.  For  example,  material  ply  failure  occurs 
at  a  point  In  the  ply  If  the  stress  state  satisfies  the  following  criterion 
[51, 

^IJ  <71  ^1  'H  ^  »  f»1  “  1*6  (1.?) 

The  Implication  of  this  type  of  criteria  Is  that  the  strength  parameters 
such  as  Fjj  and  Fj  In  (I.?)  are  constant  throughout  the  material  ply.  Conse¬ 
quently,  ply  failure  can  occur  only  as  a  stationary  event,  without  regard  to 
any  possible  modal  and/or  propagative  behaviors. 

This  Is  clearly  contradictory  to  the  observed  behavior  of  matrix  cracks 
In  laminates  as  matrix  cracking  Is  not  n  stationary  event.  Moreover,  matrix 
cracking  can  occur  as  multiple  events  vlthia  one  given  ply  at  different  ply 
stress  levels  [1,  7].  Therefore,  some  new  concept  oF  failure  Is  needel,  from 
which  a  criterion  appropriate  for  matrix  cracking  In  laminates  can  he 
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1.3  The  Concept  of  Effective  Flaws 


In  order  to  Include  the  non-uniform  ply  strength  property  In  modeling 
laminate  matrix  cracking,  one  frequently  followed  approach  Is  to  assume  a 
characteristic  strength  distribution  for  the  ply  [S],  Strength  estimation  for 
multiple  failures  In  the  ply  can  then  be  formulated  based  on  some  nrobabllls- 
tlc  stress  redistribution  arguments.  however,  this  approach  still  does  not 
relate  the  modal  or  the  propagative  nature  of  matrix  cracking. 

As  discussed  before,  the  physical  source  of  multiple  iiHtrlx  cracks  Is 
attributed  to  Inherent  material  flaws  which  distribute  randomly  throughout  the 
material  piles.  Through  the  homogeneity  assumption  made  In  ply  elasticity, 
the  physical  identity  of  the  flaws  are  lost  and  their  effects  on  ply  failure 
are  not  properly  retained. 

An  alternative  approach,  however,  is  to  represent  the  physical  material 
flaws  by  a  hypothetical  effective  flaw  distribution  at  the  macroscopic  level 
[9  -  13].  In  this  way,  the  baste  stress  analysis  scheme  can  still  he  conduc¬ 
ted  based  on  ply  elasticity,  but  the  elastic  'ledl  im  now  contains  Imaginary 
flaws. 

’^rom  this  concept.  It  Is  assumed  that  a  characteristic  distribution  of 
effective  flaws  'lavlng  random  slses  and  locations  can  be  determined  as  a 
macroscopic  property  of  the  basic  material  ply.  It  Is  required  that  the 
assumed  effective  flaws  and  their  distributional  character  must  yield  an 
effect  on  ply  failure  stmllnr  to  that  observed  In  laminate  matrl< 
cracking.  In  this  manner,  failures  la  the  Individual  piles  are  now  governe<l 
by  the  conditions  under  which  the  effective  flaws  Individually  become  nntrl< 
cracks;  the  ensuing  modal  behavior  of  crack  propagation  and  the  mul t Inl icl ty 
of  the  cracking  events  then  constlt'ite  the  observed  phenomena  at  the  nicro- 
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As  It  will  be  demonstrated  later  in  some  example  problems,  the  effective 
flaw  distribution  In  a  given  material  system  can  be  rationally  characterized 
by  a  simple  test.  The  characterization  Includes  the  form,  size  and  location 
of  each  of  the  flaws. 

When  the  form,  size  and  location  of  an  effective  flaw  In  a  laminating  ply 
are  known,  the  loathematlcal  conditions  governing  flaw  Initiation  and  propaga¬ 
tion  as  a  matrix  crack  are  then  formulated  on  the  basis  of  classical  fracture 
mechanics. 

Let  an  effective  flaw  of  known  form  in  a  considered  ply  be  of  size  a  and 
location  X.  The  stress  state  surrounding  the  flaw  can  be  determined  In  terms 
of  the  applied  laminate  stress  a  by  means  of  the  ply  elasticity  theory.  The 
critical  Incidence  at  which  the  flaw  becomes  a  matrix  crack  Is  now  derived 
from  classical  fracture  mechanics.  It  will  be  assumed  that  the  flaw/crack 
transition  Is  governed  by  the  fracture  quantity,  the  strain  energy  release 
rats  G,  at  the  flaw  tip.  Given  the  form  and  location  of  the  flaw,  G  can  he 
calculated  In  terms  of  the  applied  lamlnats  stress  o  and  the  flaw  size  a  [lA]. 
The  condition  for  the  flaw  to  propagacs  Into  a  mtrlx  crack  Is  given  by  the 
Griffith  criterion  [11]: 

G(a,  a)  =  (l.l) 

where  G^  Is  an  appropriately  detsrnlned  material  toug'iuess  for  the  particular 
node  of  the  crack  propagation. 

The  critical  lamlnats  stress  at  the  Incidence  of  flaw/crack  transi¬ 
tion  Is  determined  from  "quatlon  (l.l),  given  i  and  Gf,.  \s  for  the  prioiga- 
tlon  stahlllty,  It  Is  then  gover  ie  l  hy  the  Inequalltv  (see  also  fl'ill; 

iL  -  0 

•).a~ 
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Since  the  size  and  location  of  an  effective  flaw  are  actually  random 
variables,  and  since  the  use  of  Equation  (1.3)  In  matrlK  cracking  nrobleas  la 
limited  to  a  known  crack  Individually,  use  of  Equation  (1.3)  In  multiple 
matrix  cracks  can  lead  to  some  time-dependent  stochastic  processes.  But, 
before  a  stochastic  simulation  procedure  Is  developed,  some  basic  questions 
regarding  the  application  of  fracture  mechanics  need  to  be  discussed. 

1.4  Application  of  Fracture  Mechanics  to  Matrix  Cracks 

It  Is  well  known  that  the  classical  fracture  criteria  expressed  in 
Equations  (1.3)  and  (1.4)  apply  only  to  brittle  cracks  undergoing  stable  and 
self-similar  propagation  [13].  Whether  or  not  the  general  class  of  matrix 
cracks  considered  here  fulfills  these  requirements  needs  careful  analysis. 
Most  structural  composites  are  made  of  brittle  matrix  material,  such  as  epoxy 
resin.  Cracks  in  these  materials  can  be  considered  brittle  at  iwth  the  micro¬ 
scopic  and  the  macroscopic  scales  [111.  In  addition,  an  In-ply  matrix  crack 
will  always  remain  parallel  to  the  fibers  of  that  ply,  while  an  Iiterply 
matrix  crack,  or  a  delamlnatlon  crack,  remains  parallel  to  the  oly  Interface 
and  propagates  more  or  less  In  the  same  plane  [11].  Thus,  It  would  seem  that 
the  basic  requirements  In  applying  Equations  (1.3)  and  (1.4)  to  brittle,  self- 
similar  matrix  cracking  can  be  fulfilled  In  the  laaln. 

However,  use  of  the  fracture  criteria,  oarticularly  Equation  (1.3), 
requires  consistency  between  the  mathematical  quantity  G(o,  al  and  the 
physical  quantity  G^.  Some  details  relating  to  this  question  will  he  dis¬ 
cussed  la  the  following. 

Consider  first  the  right-hand  side  of  (1.3).  Is  a  material  ton 'hness 

property  against  matrix  rupture,  which  Is  letermlaed  by  some  suitable  phenom¬ 
enological  experiments.  \s  such,  the  >'v?jaured  far  a  given  material  system 
may  vary  wltli  a  number  of  factors.  GpecI  f  leal  Iv,  Gj.  for  unMlrecttonnl 
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fibrous  composites  changes  with  the  local  matrix-fiber  packing  geometry  as 
well  as  the  presence  of  microflaws  along  the  path  of  the  crack  [16],  But 
these  Influencing  factors  exist  at  the  microscopic  level  and  are  not  Identifi¬ 
able  at  the  macroscopic  level.  As  a  result^  the  macroscoplcally  measured 
exhibits  two  peculiar  characteristics.  The  first  Is  the  dependence  of  Gj.  on 
the  direction  of  matrix  crack  propagation  relative  to  the  ply  fiber  direction 
[11];  and  the  second  is  the  dependence  of  Gj.  on  crack  extension  Increment,  as 
the  crack  propagates  intermittently  rather  than  continually  In  real  situations 
[16].  Apparently,  both  characteristics  stem  from  the  microscopically  hetero¬ 
geneous  nature  of  the  composite  system. 

The  dependence  of  G^  on  crack  growth  direction  relative  to  the  ^Ibers  has 
been  demonstrated  In  two  experiments  of  mode-I  delamlnatlon  In  unMlrectlonal 
laminates  [17,  18].  Figure  1.6  shows  a  delamlnatlon  crack  propagating  In  the 
fiber  direction,  while  Figure  1.7  shows  a  delamlnatlon  crack  propagating 
normal  to  the  fiber  direction.  Although  the  two  types  of  delamlnatlon  crack 
occur  essentially  In  the  same  Interface  plane,  their  corresponding  fracture 
surfaces  show  different  morphology  when  viewed  at  the  microscopic 
scale.  Hence,  when  G^.  Is  deduced  from  a  macroscopic  measurement,  the  differ¬ 
ence  Is  also  considerable.  The  multl-valned  nature  of  Gj.  for  a  crack  propaga¬ 
ting  In  the  same  plane  but  In  different  directions  Is  not  generally  found  In 
conventional  metallic  materials. 

As  for  the  Intermittent,  stnp-and-go  behavior  of  'matrix  crack  propaga¬ 
tion,  this  Is  due  to  the  fact  that  the  crack  path  Is  often  wandering  around 
many  fiber  diameters.  Others  laave  also  attributed  the  phenomenon  tT  the 
viscoelastic  behavior  of  the  matrix  material.  Tn  any  e-vent,  G^  wouH  have  a 
stable  value  only  If  It  Is  averaged  over  a  finite  crack  extension  Aa.  The 
latter  can  he  as  large  as  several  fiber  diameters  and  cma  he  Interpretated  as 
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the  effective  creek  eKtenalon.  Such  e  character  Is  Inherent  to  the  fibrous 
material  system;  no  theoretical  value  for  Aa  Is  available.  Perhaps,  Aa 
discussed  here  and  AV  discussed  In  Section  1.2  have  some  common  relationship. 
Now,  the  left-hand  side  of  (1.3)  pertains  to  the  loathematlcally  computed 
G(a,  a).  The  origin  of  G  Is  traced  to  the  theory  of  Griffith  [Ih],  who  postu¬ 
lated  that  In  a  strained  elastic  body  containing  a  small  crack  of  sl^e  a,  the 
stored  elastic  strain  energy  U  will  decrease  If  the  crack  extends  stably.  The 
decreased  amount  of  the  strain  energy  Is  equal  to  the  Increased  free  surface 
energy  of  the  extended  crack.  Thus,  If  S  denotes  the  total  surface  energy  of 
the  crack,  a  balince  of  energy  at  the  Incipient  crack  extension  leads  to  the 
following  condition. 


^  ^  as 

aa  aa 


(1.3) 


In  common  practice,  the  right-hand  side  of  (1.3)  is  Interpreted  as  the 
Irrecoverable  energy  dissipated  at  the  Incipient  crack  extension  fl3l.  The 
term  would  then  Include  such  non-elastic  quantities  as  local  plastic  yielding, 
viscoelastic  dissipation,  etc.  Then,  the  resultant  luantlty  ’oecomes  neasur- 
able  only  by  phenomenological  experiment.  This  quantity  Is  comionlv  known  as 
the  critical  energy  release  rate  G^.,  which  Is  mentlond  earlier  In  F,quatlon 
(1.3). 


The  left-h’-Ki  side  of  (1.3)  Is  calculated  i.ia thematically  by  means  of  some 
elasticity  boundary-value  solution  tnethods  (e.g.  IlA]),  which  gUe  rise  to  the 
definition  of 


G(o,  a)  —  — D'l/Sa 


1  .bl 


The  mathematical  Implication  of  (1.3)  or  (l.b)  Is  than  'f  is  a  conn’nuous 
function  of  the  continuous  vsrlahle  a.  And,  at  beast  the  first  and  tie  second 
derivatives  of  the  function  G  exist  for  i  ^  0  f?ol. 
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But,  this  Implication  contradicts  the  actual  physics  of  matrix  cracks. 
As  mentioned  before,  matrix  cracks  extend  Intermittently  with  Increments  Aa 
usually  In  the  order  of  several  fiber  diameters  [161.  Furthermore,  Gj.  could 
be  measured  only  aa  an  average  over  the  finite  crack  extension  Aa. 

In  view  of  the  above.  It  seems  that  'l(o,  a)  In  (1.6)  may  be  more  appro¬ 
priate  If  defined  as 

G(o,  a)  -  _  U(aJ_Aa)_2_U^  7^ 

Aa 

and  the  stability  criteria  In  (l.A)  as 

G(o,  a  +  Aa)  -  G(o,  a)  ^  q  (1.*^) 

Aa 

where  Aa  Is  chosen  appropriately  given  the  actual  fibrous  system. 

It  Is  noted  that  between  the  definitions  in  (1.6)  and  (1.7),  the  original 
physical  meaning  Is  more  or  less  preserved.  But,  mathematically,  the  right- 
hand  side  of  (1.6)  Is  the  tangent  of  -U(a)  at  a,  while  that  In  (1.7)  Is  now 
the  secant  of  -U(a)  from  a  to  a  +  Aa;  see  Figure  l.B.  The  numerical  differ¬ 
ence  between  the  two  coulii  be  significant  depending  on  the  Form  of  'J(a)  and  on 
the  size  of  Aa. 

Clearly,  use  of  (1.7)  rather  than  (1.6)  for  the  deFlnltlon  of  G(n,  a)  can 
result  In  considerable  controversy.  For  Instance,  the  traditional  concepts  of 
crack-tip  singularity  and  stress  Intensity  Factor  K  as  a  fracture  gaantlty 
(21,  22]  may  now  require  reinterpretation. 

It  Is  known  that  the  development  )f  the  f'-conoent  was  based  on  the 
equivalence  between  y.  and  G  for  brittle  crack  In  an  elastic  contlnnn”'  121, 
22].  The  iteanlag  of  K  represtjnt®  the  latenslty  of  the  slagnlar  stress  Field 
near  the  closest  /Iclalty  of  the  crack-tip  thrt)n":h  the  general  rel  a  1 1  inshl  o 
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[211; 


aij(r,  O)  “  flj(0)  ;  l,j  =  1,  2,  3 


(1.9) 


where  Kj,  K2  and  K3  are  associated  with  the  three  modes  of  crack  propa;»atloa, 
while  (r,  0)  are  the  polar  coordinates  measured  from  the  crack-tip;  see  Fl;;ura 


The  stress  Intensity  factors  Kj  can  be  determined  In  explicit  terms  of 
the  applied  load  a  and  the  crack  size  a.  In  the  general  form 


Ki(a,  a)  =  a  gi(a) 


(I.IO) 


The  form  of  gi(a)  depends  only  on  the  geometry  of  the  crack  and  the  physical 
details  of  the  strained  body  In  consideration  (see  e.g.  (231). 

Under  the  assumption  of  Infinitesimal  local  response  at  the  crack-tip, 
the  stress  intensity  factors  Kj  are  then  related  to  the  energy  release  rate 
In  the  general  form  [20], 


G(o,  a)  =  J.jin  Nm(«»  ^ 


(l.ll) 


where  are  functions  of  the  'tiaterlal  stiffness  constants  only.  Hence, 
If  G(.  Is  a  material  property,  the  quadratic  conMoatljn  of  K^'s  must  also  be  a 
material  property. 

Use  of  K  as  a  quantity  for  crack  growth  implies  that  material  separation 
at  the  crack-tlp  region  Is  determined  solely  by  the  elastic  stress  field  which 
exists  locally;  and,  this  stress  field  Is  adequately  represented  by  the  singu¬ 
lar  function  In  the  general  form  of  (1.9). 

Now,  when  applied  to  matrix  cracking,  use  of  K  may  need  reinterpretation. 
The  reasons  are  as  follows.  First,  since  the  physical  (ustlf Icatlon  for  the 
K-concept  stems  from  the  G  -  K  equivalence,  their  mathematical  rel i t lonsh I i  as 
expressed  In  the  general  form  In  (l.ll)  becomes  tnv/alld  If  G(t,  a^  takes  on 
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the  fieflnltlon  of  (1.7)  Instead  of  (1.6).  Secondly,  since  K  sl,!»nlfles  the 
Intensity  of  a  singular  stress  field  near  the  considered  crack  tip  region, 
which  Is  often  very  small  In  slae,  say  about  a  single  fiber  llamater,  the 
meaning  of  '<  may  contradict  the  basic  requirements  of  a  finite  crack  eKtenslon 
Aa  as  well  as  the  statistically  homogeneous  postulation.  Thus,  In  *-a 

describe  the  macroscopic  behavior  of  the  composite  system,  use  of  the  K~ 
concept  for  matrix  cracking  may  also  require  some  new  Interpretation  from  both 
a  physical  and  a  mathematical  consideration. 

1.5  Cracks  In  Rl-Matarlals  Interface 

A  related  question  concerns  a  crack  that  occurs  In  the  Interface  between 
two  piles  of  different  fiber  orientations.  According  to  ply  elasticity,  this 
Is  a  crack  propagating  In  the  Interface  of  bonded  dissimilar  materials.  The 
question  relatss  to  the  stress  singularity  at  the  tip  of  the  crack,  Tn  this 
case,  the  singular  stresses  will  contain  oscillatory  components  (see 
Consequently,  the  stress  Intensity  factor  '<  or  the  strali  energy  release  rate 
G  cannot  be  computed  exactly  In  the  classical  manner.  The  oscillatory  nature 
of  the  crack-tip  stress  field  Implies  material  penetration,  which  Is  physic¬ 
ally  Inadmissible.  A  modification  of  the  classical  definition  for  K  Is 
possible  If  a  partial  crack-tip  closure  Is  Included  In  the  calculation  [241. 
Since  the  oscillatory  part  of  the  stress  field  exists  only  on  n  small  scale, 
the  Importance  of  Its  Influence  on  the  physical  quantity  G  near  the  craok-tlp 
la  usually  negligible  (20]. 

The  existence  of  crack-tip  oscillatory  stresses  Is  nttrlhut?d  to  the 
Impossibility  of  satisfying  the  r,.2ro-l.oa1  bounlnrv  conditions  near  the  ■:rack- 
tlp  point  [20] .  The  problem  seems  to  be  Inherent  to  the  theiory  of  linear 
elasticity.  The  oscillatory  part  of  the  sligular  stress  field  near  >  bl- 
materlal  crack  tip  would  disappear  IF  a  non-linear  theor/  of  elastic! tv  I^ 
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assumed  foe  the  materials  [25]. 

In  order  to  bring  some  further  clarification  to  this  problem,  a  recent 
paper  [26]  has  examined  a  case  where  an  3dhesl>;e  layer  of  finite  thickness  Is 
Introduced  In  the  Interface  of  two  dissimilar  plates.  \  delamlnatlon  crack  Is 
made  to  propagate  Inside  the  adhesive  layer.  If  the  adhesive  layer  Is  elastic 
and  homogeneous,  then  the  crack-tlp  stress  singularity  will  contain  no  oscill¬ 
atory  components.  The  strain  energy  release  rate  G  in  this  case  can  be  compu¬ 
ted  In  the  conventional  manner  provided  the  various  material  Interface  condi¬ 
tions  are  Incorporated  In  the  computation.  When  the  adhesive  layer  thickness 
Is  reduced  to  nearly  aero,  the  computed  G  uniformly  approaches  a  finite 
quantity,  see  ^'Igure  I.IO.  It  would  thus  seem  that  the  oscillatory  stress 
field  of  a  hl-materlal  crack  Is  not  an  Issue  for  delamlnatlon  analysis, 
particularly  In  view  of  the  requirement  of  the  statistically  homogeneous  body 
concept. 

As  demonstrated  In  [26],  for  a  delamlnatlon  crack  between  tvo  piles  of 
different  fiber  orientations,  the  mixed-mode  energy  release  rate  n  can  also  be 
adequately  calculated  by  a  finite  element  method  which  dlsregarls  the  stress 
singularity  altogether. 

1.6  The  Finite  FAement  ^r vice  lures 

From  t'ne  classical  fracture  mechanics  viewpoint,  both  the  G-concent  and 
the  K-concept  are  equally  applicable  to  brittle  crac.k  propagation  problems. 
Relative  ease  la  computation  Is  oft»n  the  reason  dictating  which  one  >F  the 
concepts  Is  employed.  In  the  case  of  crack  In  loetals,  the  "-Ciinc'ot  has 
traill tlonally  been  prefarreil. 

When  a  crack  Is  situated  In  a  body  of  complex  boundaries,  nui'-loal 
procedures  such  as  the  contemporary  finite  element  technliues  are  frequentlv 
used.  In  this  case,  the  computation  of  G  Is  usually  slmpl'.er  than  the  como'ita- 
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} 

J  tlon  of  K.  In  either  case,  uncertainties  over  numerical  accuracy  often  arise 

I  because  of  the  numerical  technique  used  to  approximate  the  slnt^ular  stress 

field. 

These  concerns,  however,  are  largely  removed  If  the  discrete  form  for  G 
I  as  expressed  by  (1.7)  Is  computed.  But  the  practical  question  remains  as  to 

the  proper  definition  for  the  effective  crack  extension,  Aa.  This  question  Is 
related  to  the  Inherent  character  of  the  physical  material  system.  It  can  be 
I  answered  only  by  experience  and  experiment.  Except  for  the  uncertainty  over 

the  choice  of  Aa,  the  finite  element  technique  is  a  useful  and  direct  means 
,  for  computing  G,  as  expressed  In  the  form  of  (1.7). 

I  Figure  l.ll  shows  the  finite  element  representation  of  a  plane  crack  In  a 

2-dlmenslonal  solution  field.  Here,  the  tip  of  the  crack  of  length  a  Is  at 
node  c.  Given  the  applied  far-fleld  load  o,  the  entire  elastic  stress  field 
I  can  be  determined  within  the  framework  of  the  finite  element  method.  The 

solutions  not  only  Include  the  stresses  and  displacements  (at  the  nodes)  but 
also  the  total  strain  energy  stored  la  the  body,  U(a).  Then,  by  a  properlv 
^  chosen  Incremental  crack  extension  Aa  la  the  direction  oF  the  original  crack, 

the  node  c  Is  opened  up  Into  two  separate  no<les  f  and  g,  whereby  the  crack  tip 
advances  to  node  h.  With  this  new  crack  configuration  and  by  keeolag  the  far- 
^  field  conditions  unchanged,  a  new  solution  Is  obtained  which  yields  the 

displacements  (u,  v,  w)  for  nodes  f,  g  and  h,  as  well  as  the  stared  strain 
energy  U(a  +  Aa) .  This  stan  readily  yields  the  strain  energy  release  rst' 
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n) 

as 
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The  work  done  In  closing  the  crack  extenslan  ^n  Is  then  given  b.- 
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AW  =  1/2  iF^Cuf  -  Ug)  +  Fy(vf  -  Vg)  +  F2(wf  -  Wg)l  (1.12) 

According  to  the  work-energy  equivalence  principle  [22],  the  strain 
energy  release  rate  during  the  finite  crack  extension  Aa  can  he  alternatively 
defined  as 

G(‘^,  a)  =>  AW/Aa  (1.11) 

Expression  (1.13)  gives  the  same  results  as  (1,7).  8ut,  the  advantage  of 
the  expression  In  (1.13)  Is  that  It  also  gives  the  modal  components  In  0  If 
the  crack  propagation  Is  In  mixed-mode.  Thus,  by  way  of  (1.12),  Equation 
(1.13)  takes  the  form, 

G(a,  a)  =  Gj(a,  a)  +  G];[(o,  a)  -t  G-[£[((;,  a)  (1.14) 

where 

G-[(a,  a)  =  F2  (wf  -  Wg)/2Aa 

Gu;(o>  a)  =  Fy  (vf  -  Vg)/2Aa  (I.IS) 

Gu;[(o,  a)  =  F^  (uf  -  Ug)/2Aa 

By  successive  opening  and  closing  of  the  nodes  along  the  crack  propaga¬ 
tion  path,  the  quantity  G(a,  a)  can  be  calculated  as  a  function  of  the  crack 
length  a. 

This  technique  was  apparently  first  applied  to  delamlnatlon  cracks  In 
laminates  by  Ryblckl  and  Kannlnen  [27].  It  was  later  extended  to  slnulite 
transverse  cracks  [9]  and  contoured  delamlnatlon  la  3-d 1 aens tonal  stress  field 
by  Wang,  et.  al  [28|. 

1.7  Summ.3  ry 

In  this  chapter,  several  basic  i)hyslcal  and  mathcaatlcal  coiacepta  frav e 
been  discussed.  These  are  pertinent  t)  a  consistent  development  of  a  failure 
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analysis  inethodology  for  niatrlx  cracks  In  composite  laminates.  At  this  point, 
a  summary  discussion  of  these  concepts  Is  briefly  presented  below: 

(1)  Matrix  cracks,  their  Initiation,  propagation  and  interaction  consti¬ 
tute  a  complex  failure  process  In  fibrous  composite  materials,  ^or  the  class 
of  laminates  that  are  made  by  stacking  unldlrectlonally  fiber  reinforced  plies 
of  different  fiber  orientations,  the  matrix  cracking  patterns  are  generic  to 

the  lamination  geometry  and  the  nature  of  the  applied  load.  The  general  | 

cracking  process  Is  seen  to  follow  a  certain  mechanistic  rule,  and  the  problem 

appears  to  be  amenable  to  some  form  of  analytical  description.  ^ 

(2)  When  viewed  at  the  macroscopic  scale,  a  matrix  crack  can  be  classl-  ) 

fled  as  one  of  two  fundamental  forms.  The  In-ply  crack  which  separates  | 

1 

through  the  thickness  of  a  ply  Is  confined  to  propagate  parallel  to  the  fibers  j 

In  that  ply.  The  Interply  crack,  known  also  as  delamlna t Ion ,  is  propagated  In  f 

a  piano  parallel  to  or  Inside  of  ply  Interfaces.  Both  mav  be  considered  as  j 

J 

self-similar  propagating  cracks  at  the  macroscale.  j 

(3)  When  viewed  at  the  microscopic  scale,  l.e.  with  precise  flbe: — matrix  j 

structure  taken  Into  account,  matrix  cracks  are  thought  to  originate  from  ' 

I 

material  flaws  which  exist  randomly  either  In  the  matrix  phase  or  on  the 

I 

« 

fiber-matrix  Interfaces.  The  Individual  propagation  and  mutual  coalescence  of  ^ 

these  flaws  lead  to  matrix  cracks  of  large  proportion,  which  are  then  observe-  ^ 

able  at  the  macroscale.  The  Initiation  and  growth  behaviors  oF  matrix  cracks 

c 

I 

at  the  macroscale  depend  Intrinsically  on  the  local  details  of  flaw  growth  and  / 

flaw  coalescence.  These  local  details,  however,  cannot  he  practically  j 

modelled  because  of  the  multiplicity  and  the  nrohahl 1  Is 1 1 c  nature  iif  the  J 

material  flaws.  J 

(A)  The  accepteil  theoretical  basis  for  stress  md  Fnllnre  anal/sls  of  j 

'i 

composite  laminates  has  been  the  concent  of  s ta 1 1 s 1 1  ml  1 v  homo Tcneous  body.  > 

I 
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From  this  coocept,  the  theory  of  ply  elasticity  has  traditionally  been  derived 
and  applied.  But,  this  theory  compromises  the  presence  and  the  effects  of  the 
material  flaws  on  failure.  Consequently,  these  microscopic  effects  cannot  be 
recovered  rationally  In  the  ensuing  stress  and  failure  analyses  In  the  conven¬ 
tional  manner. 

(5)  In  order  to  overcome  the  predicament  caused  by  micro-macro  translo¬ 
cation  of  the  material  flaws,  an  alternative  approach  Is  followed.  The 
approach  Is  to  work  within  the  realm  of  ply  elasticity,  but  to  Include  the 
presence  and  the  effects  of  the  material  microflaws  In  the  form  of  assumed 
effective  flaws.  Specifically,  given  a  finite  volume  of  the  basic  ply  'nater- 
lal.  It  assumes  that  a  charocterlstlc  distribution  of  effective  flaws  can  be 
Identified  which  has  the  equivalent  effects  as  the  Intrinsic  flaws  that  exist 
at  the  microscale.  Although  the  existence  of  the  effective  flaws  Is  hypothe¬ 
tical,  they  nevertheless  provide  a  viable  alternative  for  modelling  matrix 
cracks  at  the  inacroscale.  As  such,  the  assumed  effective  flaw  distribution 
represents  a  macroscopic  ply  property. 

(6)  The  Introduction  of  the  assumed  effective  flaw  distribution  not  only 
retains  the  probabilistic  and  the  multlpllstlc  charictarlstlcs  of  the  laatrlx 
cracking  process  as  a  whole,  It  also  enables  the  description  of  the  propaga¬ 
tive  behavior  of  the  Individual  matrix  cracks  themselves.  Specifically,  the 
Individual  flaw  Is  treated  as  a  crack  whose  propagation  mechanisms  Is  governed 
by  the  criteria  derived  from  the  theory  of  linear  elastic  fracture  mechanics; 
the  collective  behavior  of  a  distribution  of  random  flaws  Is  then  treated  by  a 
method  of  stochastic  processes. 

(7)  Within  the  logical  context  ootllneil  above,  two  problems  of  fundamen¬ 
tal  nature  arise.  The  first  Is  related  to  the  actual  characterization  of  the 
effective  flaw  distribution,  given  a  finite  volnnc?  i>f  the  plv  material.  By 


NADC-85118-60 


definition,  this  distribution  can  be  characterized  only  by  phenoiiienolo.f»lcal 
means.  And,  as  such,  a  considerable  amount  of  empiricisms  would  be  Involved 
In  the  characterization.  This  will  undoubtedly  pose  uncertainties  as  to  the 
uniqueness  and  the  reliability  of  the  determined  flaw  distribution  as  a  ply 
material  property. 

The  second  problem  concerns  the  consistency  between  the  physics  and  the 
mathematics  associated  with  the  crack  propagation  criteria  used.  On  one  hand, 
the  theory  of  elastic  fracture  inechanlcs  assumes  that  crack  propagation  is 
governed  by  conditions  which  exist  in  the  closest  vicinity  of  the  crack-tip, 
where  the  dimension  of  the  region  of  interest  Is  often  smaller  than  the 
diameter  of  a  single  fiber,  or  the  actual  size  of  a  microflaw.  On  the  other 
hand,  the  actual  physics  of  'oatrlx  cracking  seems  to  Indicate  that  nlcroscoplc 
activities  near  the  crack-tip  may  extend  over  a  region  larger  than 

several  fiber  cross-sections. 

(8)  To  alleviate  these  difficulties,  a  modification  of  the  classical 
criteria  of  crack  propagation  Is  proposed,  Namely,  only  the  effective  frac¬ 
ture  quantities  are  defined  and  used.  In  particular,  a  discrete  form,  or  an 
averaged  value  of  the  strain  energy  release  rate  G  over  a  finite  crsck  .exten¬ 
sion  Aa  Is  defined,  see  Equation  (1.7).  The  chol.ce  for  the  proper  size  of  the 
crack  extension  An  will  depend  on  such  factors  as  the  size  of  material  hetar.j- 
genelty  (a.g.  fiber  diameter,  mlcr.aflaws,  etc.)  l.a  the  actual  composite  system 
and  the  lamination  fabrication  variables. 

Using  the  concepts  oreseated  ab.')''e,  the  f^ractur.e  descrl.ptlon  far  ii.atrl< 
cracks  'iiay  be  .auj'inented  by  a  numerical  finite  el.»ii.?nt  simulation.  Agalast 
this  background,  an  aa.alytl.:al  meChoflology  Is  nreseuti'l  and  applied  ti  se-'eral 
natrlx  cracking  probl.em;;.  Concurrent  experiments  have  been  c.induo t si  to 
adjudicate  the  theory.  Most  if  tije  ditslls  are  presented  la  the  followlag 
chapters. 
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ease  Pate  as  a  Function  of  the  Thickness  of 
,  1  Interface  1761. 
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CHAPTER  II.  ANALYSIS  OF  INTRAPLY  CRACKING 

2.1  Intraply  Versus  Interply  tlatrtx  Cracks 

As  was  seen  In  the  previous  chapter,  structural  composites  made  from 
laminated  unidirectional  piles  frequently  suffer  matrix  cracks  well  before 
final  laminate  rupture.  These  matrix  cracks  can  develop  Into  dama^re  zones  of 
various  degrees  of  severity  under  service  loading. 

To  visualize  the  physical  nature  of  matrix  cracks  at  the  macroscale,  It 
Is  useful  to  examine  the  sketch  In  Figure  2.1.  Individually,  each  unidirec¬ 
tional  ply  Is  highly  anisotropic.  For  example,  the  ply  stiffness  In  the  fiber 
direction,  E^,  can  be  orders  of  magnitude  larger  than  the  stiffness  normal  to 
the  fiber  direction.  Eg-.  Similar  differences  also  exist  lietween  the  In-ply 
Poisson  ratios  and  the  thermal  expansion  coefficients  In  the  longitu¬ 

dinal  direction  qcl,  and  the  transverse  direction,  aj,  etc.  Table  II-l  shows 
the  disparity  between  the  fiber  dominated  and  the  matrix  dominated  effective 
ply  properties  of  a  typical  graphite-epoxy  unidirectional  system. 

Clearly,  the  In-plane  property  mismatch  bet-’een  plies  of  different  fiber 
orientations  can  cause  complex  in-ply  and  Interply  stresses  when  the  piles  are 
lamlnatsd  together.  The  thermal  expansion  mismatch  between  plies  In  laminates 
cured  at  high  temperature  will  result  In  considerable  thermal  residual 
stresses  upon  cooling  to  ambLent  temper  a  tire.  Likewise,  a  mismatch  of  stiff¬ 
ness  and  Poisson  ratio  between  piles  complicates  the  laminate  stress  field 
when  external  service  load  Is  applied. 

Since  the  ply  tensile  strength  transverse  to  the  fibers  Is  usuallv  orlers 
of  magnitude  smaller  than  the  strength  la  the  finer  llrectlon,  (sec  n.g.  Table 
II-l),  matrix  cracks  due  to  la-ply  tenslan  s'l'  or  literolv  tenslan  frequent¬ 
ly  occurs  first.  Rut,  whether  an  Intraply  matrix  (’.rick  occurs  h^'^ore  an 
laterply  'natrlx  crack  In  a  given  laml'aata  denends  on  the  laminnt?  geoaetrv  and 

II 
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tlie  exact  nature  of  loading.  Particularly,  geometrical  factors  profoundly 
Influence  the  laminate's  In-ply  and  Intarply  stress  fields. 

To  Illustrate,  compare  the  following  two  cases  between  the  [0/90] 3  and 
the  [+25/90] 3  laminate  coupons: 

The  [0/90)3  Tensile  Coupon.  Assume  that  the  [0/90)3  laminate  coupon  is 
long  and  Is  loaded  uniformly  In  axial  extension.  Onder  such  a  condition,  the 
laminate  stress  field  Is  Independent  of  the  axial  coordinate  x  and  the  three- 
dimensional  laminate  stresses  can  he  analyzed  by  a  two-dimensional  generalized 
plane  strain  formulation  In  the  y-z  plane.  This  problem  has  originally  been 
treated  by  Pipes  and  Pagano  [30];  and  a  later  treatment  using  a  sepclal  finite 
element  method  was  carried  out  by  Wang  and  Crossman  [31). 

Now,  let  the  material  system  be  T300/934  graphite-epoxy  system  (see  Table 
II-l  for  ply  properties)  and  consider  the  stress  field  Induced  by  an  applied 

I 

far-fleld  strain,  e^  =  I  po.  The  dominant  In-ply  and  Interply  stresses  In  the 

OO^-ply  are  sketched  In  Figure  2.2.  It  Is  seen  that  the  tensile  stress  Is 

due  to  the  applied  tension,  while  a  compresslwe  oy  stress  Is  caused  by  the  ^ 

Poisson  ratio  mismatch  between  the  •')‘’-ply  and  the  90'’-ply.  On  the  other  hand, 

the  tensile  0^  and  the  shearing  Ty^  stresses  are  effecte>1  by  and  concentrated 

near  the  laminate  free  edge  [30],  figure  2.3  shows  the  a.;  distribution  on  the  j 

laminate  mid-plane  (90/90  Interface)  as  a  function  of  distance  from  the  free 

edge.  Note  the  sharp  rise  li  magnltu<le  for  near  the  laminate  free  edge. 

In  this  example,  the  numerical  value  of  the  la-ply  stress  la  the  90°- 
ply  Is  about  11.7  Kpa  per  ije,  while  the  maxlmun  i;.  stress  near  the  free  edge 
Is  less  than  2  Kpa  per  pe.  It  Is  clear  by  comparison  that  In-ply  matrix 
cracks  (transverse  cracks)  cmised  by  ire  likely  t.>  occur  la  the  "•O^’-ply  ' 


before  the  ed.ge  stress  Induced  lataroly  delaml nation. 
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The  [-25/90]a  Tensile  Coupon.  Now,  consider  the  [rtS/golg  coupon  under 
similar  situation.  In  this  case,  the  +25“-plles  as  a  unit  has  a  smaller  axial 
stiffness,  but  has  a  larger  In-plane  Poisson  ratio  v^y.  Consequently,  both 

the  In-ply  and  the  Interply  stresses  near  the  free  edge  region  are  consider¬ 

ably  altered,  especially  the  Interply  stress  oz  near  the  free  edge  region 
which  Increases  dramatically  In  value.  As  seen  In  Figure  2.4,  the  value  of  oz 
at  the  free  edge  Is  about  20  Kpa  per  ye,  a  ten-fold  Increase  from  the  previous 
example. 

On  the  other  hand,  the  la-ply  tensile  stress  In  the  OO^-ply  remains 
about  the  same  at  11.7  Kpa  per  pe.  Hence,  It  may  be  Inferred  that  Interply 

cracks  In  the  form  of  edge  de lamination  In  the  90° /90®  Interface  occurs  before 

the  transverse  cracking  In  the  90°-ply. 

As  both  the  In-ply  matrix  crack  and  the  Interply  matrix  crack  ace  local¬ 
ized  failures,  the  occur  ranee  of  one  may  precipitate  the  occurrence  of  the 
other.  In  the  following,  however,  the  basic  mechanisms  of  these  two  basic 
forms  of  matrix  cracking  s'lall  be  examined  separately. 

2 . 2  Physical  Characteristics  of  Intraply  Cracking 

Consider  first  the  physical  charseters  of  Intraply  matrix  cracking. 
Take,  for  example,  the  [0/90] 3  type  tensile  coupon  discussed  above,  because 
of  the  relatively  small  Interlaminar  edge  stresses,  formation  of  matrix  cracks 
Is  largely,  if  not  exclusively,  inside  the  thickness  of  the  90°-plles  la  the 
form  of  transverse  cracks.  It  Is  thus  possible  to  use  this  example  and 
examine  the  cracking  events  without  regacl  to  the  influence  of  laminate  free 
edge  effects. 

The  Onset  and  Growth  Char ic ter ts tics.  When  viewed  at  the  macroscopic 
scale,  the  formation  of  a  transverse  crack  In  the  90°-pl/  Is  simply  a  sudden 
separation  of  the  collimated  fibers  that  lie  normal  to  the  applied  tension. 
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Within  a  single  specimen,  numerous  such  cracks  could  form  as  a  sequence  of 
load-dependent  events;  onset  of  the  first  crack  occurs  when  the  applied  load 
reaches  a  certain  critical  value,  and  similar  successive  cracks  emerge  as  the 
applied  load  Increases.  Experimental  records  by  x-radlographs  are  shown  la 
Figure  2.5  for  a  [0/90] g  coupon  under  monotonlc  tension.  The  pictures  show 
the  emergence  of  transverse  cracks  as  a  function  of  the  applied  load.  Erom 
the  photographs  of  this  kind,  a  plot  of  the  transverse  crack  density  (cracks 
per  unit  coupon  length)  versus  the  applied  tension  can  be  obtained. 

Figure  2.6  shows  a  family  of  such  experimental  plots  from  tests  of 
T300/934  [0/90ji/0]  laminate  series  with  n  =  I,  2,  3  and  4  [32].  Except  for 
the  laminate  of  n  =  1  for  which  no  visible  transverse  crack  was  detected 
before  final  laminate  rupture,  all  three  cases  (n  =  2,  3,  4)  yielded  distinct¬ 
ive  transverse  crack  density  vs.  applied  tension  relationships.  This  data 
shows  that  the  90°-ply  transverse  cracking  process  possesses  both  a  probabil¬ 
istic  and  a  deterministic  character.  Specifically,  the  occurrence  of  any  one 
individual  crack  Is  highly  probabilistic,  while  the  overall  crack  development 
process.  Including  the  Initiation  load  and  the  crack  density  vs.  load  rela¬ 
tionship  as  a  whole  are  deterministic. 

Several  geometrical  factors  were  found  to  Influence  the  characteristics 
of  the  cracking  process  [32],  In  addition  to  the  limlnate  stacking  sequence, 
the  actual  thickness  of  the  bO^-plles  (or  the  number  n)  also  Influences  the 
cracking  behavior.  This  can  be  seen  from  the  data  la  Figure  2.6,  whldi  shows 
that  a  thinner  90°-layer  Is  capable  of  reaching  a  higher  crack  density  than  a 
thicker  90°-layer.  Also,  the  critical  tensile  strati  Cf,,.  for  the  onset  )f  the 
transverse  crack  process  decreases  with  Increased  bO^-layer  thickness,  see 


Table  Tl-2. 
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The  Kinematics  oC  Crack  Propagation.  In  order  to  develop  a  descriptive 
model  for  Intraply  cracking,  the  kinematics  of  the  Individual  crack  propaga¬ 
tion  must  first  be  understood.  Despite  numerous  experiments,  the  exact  nature 
of  the  00®-ply  transverse  crack  propagation  is  unclear.  From  inacroscoplc 
observations.  It  has  been  suggested  that  the  crack  might  Initiate  from  a 
cluster  of  flaws  which  lie  thlckness-wlse  In  the  DQ^-ply.  As  load  reaches  the 
critical  value,  a  sudden  coalescence  occurs  accompanied  by  audible  acoustic 
emissions  (see  e.g.  (33]).  This  dynamic  character  may  be  due  to  the  very 
rigid  axial  stiffness  of  the  O^-plles  constraining  the  flaws  In  the  hO^-ply  to 
propagate  primarily  In  the  thickness  direction.  The  propagation  Is,  however. 
Immediately  arrested  or  blunted  at  the  0/90  Interface,  see  Illustration  In 
Figure  2.7. 

The  thlckness-wlse  crack  propagation  is  merely  con]ecture,  though  it 
conforms  to  the  phenomenonologloal  behavior  observed  experimentally.  This 
conjecture  provides  a  basis  for  defining  and  characterizing  the  "effective" 
flaws,  as  well  as  for  simulating  the  multiple  cracks  propagation  process. 

2 . 3  Representation  of  the  Effective  Flaws 

For  practical  purpose,  assume  that  the  In-ply  "effective"  flaws  are  one- 
dlmenslonal,  being  orientated  In  the  ply  thickness  direction  as  Illustrated  In 
Figure  2.8.  The  linear  size  of  an  Individual  flaw  Is  denoted  by  2a,  which  Is 
bounded  by  the  total  thickness  of  the  90'’-plles,  2nt.  The  relative  spacing 
between  any  two  adjacent  flaws  Is  denoted  by  a,  which  can  he  detecnlael  by 
their  relative  positions  (xjfj  -  xj).  leace,  both  a  and  s  are  real,  nonltlve 
random  variables.  liven  the  laminate,  the  distribution  Function  f(a)  inl  the 
distribution  function  f(s)  are  sufficient  to  define  the  charscterl sties  of  the 
Flaws  In  the  90°-plles. 

As  will  be  shown  next,  the  distribution  function  f ( a )  and  f(s)  in 
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laminates  with  grouped  90®-plies  can  be  obtained  in  terms  of  the  flaw  distri¬ 
bution  functions  of  a  single  90®-pl7  through  a  statistical  v; olume- theory . 
Then,  It  Is  necessary  only  to  characterize  the  flaw  distribution  functions  of 
the  single  90“-ply  as  a  basic  ply  property. 

Effective  Flaws  In  Single  Material  Ply.  Let  2t  be  the  nominal  thickness 
of  the  basic  90®-ply.  For  a  unit  length  of  the  90“-ply,  say  1  cm  long,  let 
there  be  M  effective  In-ply  flaws  which  lie  thlckness-wlse  and  distribute 
randomly  along  the  1  cm  length.  Figure  2.8.  If  the  half  size  of  the  Ith  flaw 
Is  denoted  by  (ao)i  end  Its  location  Is  denoted  by  (xq)i,  then  the  sets 

I  “  1»  h  characterize  the  In-ply  flaws  of  the  90®-ply  of 

unit  length  (1  cm). 

Now,  let  the  set  {(ao)i}  be  represented  by  a  Welbull  distribution  of  the 
continuous  variable  a^  In  the  form  [34], 

Fi(ao)  =*  1  -  exp  (-(ao/e^)**®!  (2.1) 


where  and  Pq  are  the  shape  and  the  scale  parameter  of  the  cumulative 
distribution  Fi(ao). 

In  the  mathematical  form  of  (2.1),  the  variable  a^  extends  from  0  to 
But,  as  a  practical  matter,  Sg  must  be  bounded  between  0  and  t.  Por  Instance, 
the  l-percentlle  of  Sg  should  be  greater  than  0  and  the  99-percentlle  of  ag  be 
smaller  than  t.  Then  the  following  constraints  on  the  values  of  og  and  Pg  may 
be  Imposed: 

l/a 

[(0.01)  '’1  Pg  >  0 


[(4.6) 


<  c 


(2.2) 


The  charncterlzatlon  of  the  In-ply  flaw  sl»:e  distribution  then  rsducss  to 
determining  the  number  M  and  the  parsmetirs  xg  and  pg,  subjected  tn  the 
possible  constraints  In  (2.2). 
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As  It  Mill  be  showa  next,  the  distribution  functions  f(a)  and  f(s)  In 
laminates  with  grouped  90“-plles  can  be  obtained  In  terms  of  the  flaw  dlstrt- 

As  for  the  flaw  spacing  distribution  f(8)  or  the  location  sat  {(xqH)* 
I  =  1,  M,  It  Is  often  sufficient  to  select  a  suitable  random  number  generator 
that  gives  a  set  of  M  random  valties  between  the  Interval  (0,  1).  These  random 
values  then  constitute 

Effective  Flaws  In  Grouped  Piles.  In  laminates  of  more  than  one  <i0°-ply 
grouped  together  as  In  the  [0/90^/0!  family  when  n  >  2,  the  effective  flaws  In 
the  grouped  90‘’-plles  are  not  the  same  3l^e  as  those  In  the  single  'i0‘’-pl/,  \ 

simple  scheme  loay  be  used  to  determine  the  effective  flaw  sizes  In  the  grouped 
90“-plles,  given  the  flaw  size  distribution  of  the  single  90“-ply. 

The  scheme  la  based  on  a  statistical  consideration  of  Increased  volume 
[34] .  It  may  be  perceived  at  the  first  glance  that  the  grouped  90“-plles 
would  behave  as  n  Individual  plies  arranged  In  parallel,  similar  to  a  bundle 
of  threads  loaded  In  uniform  tension  [351.  However,  because  the  piles  ara 
laminated  together,  failure  of  any  one  90°-ply  may  lead  to  unstable  propaga¬ 
tion  through  the  thickness  of  the  grouped  90'’-plles.  Consequently ,  failure  of 
the  grouped  90'’-plles  more  likely  behaves  as  n  piles  arranged  In  series. 

Assume  that  for  every  unit  length  (I  cm  long)  of  the  hO^-ply,  there 
exists  a  discrete  tensile  strength  distribution,  {(ao)i}  I  =  1,  M;  the  set 
{(oo)l}  l3  assumed  to  be  related  one-to-one  with  the  effective  flaws  {(aQ)j} 
through  the  linear  fracture  ii»echantcs  law: 

(oo)l  =  K//Ta^  ;  I  =  1,  M  (2.3) 

where  K  Is  some  constant. 

Now,  let  the  {(oo)i}  set  he  represented  by  a  Welbull  failure  probability 
function  of  the  continuous  variable  oq, 

Pl(ao)  -  1  -  exp  [-(no/, go)  o] 

37 


'  7 


.  k,^  «•»» 
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The  parameters  oq  Po  *■**  (2.4)  are  Implicitly  related  to  oq*  8o 
In  (2.1)  through  the  relationship  of  (2.3), 

When  n  90“-plle8  grouped  together,  the  probability  of  failure  based  on  a 
In-serles  arrangement  la  given  by, 


Pn(ao^  “  1  -  (1  -  Pi(oo)l" 


(2.S) 


which,  upon  substitution  of  Pi(ao)  In  tbe  form  of  (2.4),  leads  to  the  Welbull 


function: 


where 


Pn(<Jo^  =  1  -  exp  [-(co/Bn^^*"! 


On  ”  ^  ?a  =■  Sod/n) 


(2.6) 


(2.7) 


The  second  result  in  (2.7)  suggests  that  the  charecterlstlc  failure 
parameter  of  n  90”-plles  in  group  is  determined  by  ^  and  p©  of  the  single 
ply  strength  distribution  in  (2.4).  In  fact,  the  same  result  for  can  be 
obtained  If  the  grouped  plies  are  assumed  to  behave  like  a  bundle  arranged  in- 
parallel  if  the  number  n  In  the  bundle  Is  small,  say  n  <  6  (36). 

In  view  of  the  specific  relationship  between  the  strength  oq  ^nd  flaw 
size  a^.  Equation  (2.3),  It  may  be  assumed  that  the  effective  flaw  size 
distribution  In  grouped  plies  obeys  a  volume-rule  similar  to  (2.7).  Thus,  let 
the  effective  flaws  of  n  hO^-plles  in  group  be  represented  by  the  random 
variable  1^,;  and  a^  is  determined  from  and  tjQ  of  the  single  ply  by  the 


relationship. 


2/a 

an  =•  ao(n) 


(2.8) 


As  for  the  location  distribution  {(’Cj,)i},  It  Is  assumed  unchanged  for  any 
number  of  90'’-plles  grouped  together. 


b"  V  •-*  V* 
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It  should  be  emphasized  that  the  volume-law  (2.8)  for  effective  flaw  size 
distribution  In  grouped  90"-plles  Is  only  empirical  and  should  therefore  be 
regarded  as  tentative.  Its  usefulness  remains  to  be  experimentally  verified. 
Further,  research  based  on  microscopic  examination  and  statistical  mechanics 
Is  needed  In  order  to  define  a  flaw  representation  scheme  on  more  rigorous 
grounds. 

2.4  Shear-Lag  Effect  on  Crack  Spacing 

When  a  90“-ply  transverse  crack  Is  formed,  the  local  tensile  load 
formally  carried  by  the  QO^-ply  Is  now  transferred  to  the  adjacent  O^-plles. 
Because  of  the  Interface  bonding  strength  between  the  0“  and  the  00“  piles,  an 
Interlaminar  shear  stress  Txz  Is  developed  near  the  transverse  crack  termini, 
see  Figure  2.9.  This  shear  stress  Is  In  fact  singular  at  the  crack  root  but 
decays  exponentially  a  distance  away.  As  the  shear  stress  decays,  the  In-sltu 
tensile  stress  ox  lo  the  90“-ply  regains  Its  far-fleld  magnitude.  The  ’’size" 
of  this  local  load-transfer  zone,  known  also  as  the  shear-lag  zone.  Is  roughly 
proportional  to  the  total  thickness  of  the  O0“-plles  which  contain  the  trans¬ 
verse  crack.  For  the  [0/901 g  type  tensile  coupons  considered  In  this  section, 
the  shear-lag  zone  Is  about  8  times  the  total  thickness  of  the  90“-plles 
[37]. 

The  effect  of  the  90“-ply  transverse  crack  Is  that  It  alters  the  stress 
state  locally.  If  a  90“-ply  flaw  Is  located  near  the  transverse  crack,  or  lie 
Inside  Its  shear-lag  zone,  the  flaw  will  become  less  likely  to  Initiate  and 
propagate  Into  a  crack.  If  the  flaw  Is  located  outside  the  shear-lag  zone, 
then  the  presence  of  the  neighboring  transverse  crack  will  not  have  any  effect 
on  the  behavior  of  the  flaw.  Kence,  the  shear-lag  effect  Indirectly  regulates 
the  spacing  of  transverse  cracks.  This  effect  can  he  seen  from  the  crack 
density  versus  load  relationships,  e.g.  shown  In  Figure  2.6. 
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Clearly,  It  ts  Important  to  Include  this  shear-lag  effect  either  Impli¬ 
citly  or  explicitly  In  the  predictive  model  for  multiple  cracks. 

2.5  Formulation  of  the  Simulation  Model 

The  essence  of  the  predictive  model  Is  that  the  transverse  cracking 
process  is  regarded  as  a  sequence  of  local  events  which  occur  during  the 
course  of  loading.  Each  event  la  a  fracture  propagation  problem  Involving  the 
growth  of  a  certain  effective  flaw  Into  a  transverse  crack.  The  conditions 
governing  the  Individual  flaw  growth  depend  on  the  local  stress  state  surroun¬ 
ding  the  flaw.  This  also  Includes  the  effects  of  shear-lag  zones  and  the 
concurrent  cracking  state  In  the  90“-ply.  Since  the  sizes  {ai}  and  locations 

{x^}  of  the  flaws  are  random,  the  exact  order  of  the  cracking  events  Is 

actually  stochastic  In  nature. 

In  what  follows,  a  probabilistic  simulation  procedure  derived  from  the 
Monts-Carlo  search  technique  Is  presented  for  the  transverse  cracking  process. 

Onset  of  the  First  Crack.  Consider  the  [')/90]3  type  laminate  coupon 
under  uniaxial  tension.  \ssume  that  there  ar?  M  effective  flaws  In  a  unit 
length  of  the  grouped  ^O^-plles.  The  sizes  of  the  Flaws  (=  23^)  are  represen¬ 
ted  by  the  sat  {ai}  and  their  locations  by  {xj};  I  =  1,  M.  In  this  case,  all 

flaws  are  one-dlmenslonal  and  lie  thlckness-wlse  In  the  9r)®-plles,  as  shown  la 
Figure  2.8. 

Under  the  applied  laminate  tension,  say  a  uniform  tensile  strain  the 

largest  flaw  in  {aj},  a^ax*  ’•'HI  the  first  to  propagate  Into  a  transverse 
crack  located  at  <*.  The  critical  value  of  (^x^cr  this  to  happen  thus 
defines  the  onset  load  of  the  transverse  crack  events. 

The  determination  of  (ex)cr  essentially  a  Fracture  mechanics  problem. 
As  postulated  earlier,  the  Individual  effective  flaw  would  behave  as  a  small 
crack.  Under  the  applied  tension  ex  and  possibly  also  a  uniform  tenperatire 
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drop,  at,  due  to  Laminate  curing,  the  flaw  of  size  2ai  would  undergo  mode-I 

o 

propagation.  The  energy  release  rata  which  Is  available  at  the  flaw-tip 
can  be  calculated  by  means  of  the  2-dlmen8lonal  generalized  plane  strain 
finite  element  crack-closure  procedure  as  outlined  In  [11].  Thus, 

o  _  _  _  _  . 

(ex,AT,ai)  =  (/Cg  •  e^  +  /Ct  •  AT]2  (2t)  (2.9) 


where  Cg,  and  are  coefficients  depending  on  the  value  of  a^  and  the  proper¬ 
ties  of  the  particular  lamlnata  In  question;  they  are  Independent  of  the  loads 
ejj  and  aT.  The  quantity  (2t),  which  Is  the  thickness  of  one  ply,  represents 
the  linear  scale  of  the  finite  element  model.  That  Is,  In  the  finite  element 
model,  ply  thickness  Is  set  to  a  unity. 

Onset  of  flaw  propagation  l.s  governed  by  the  criterion  (1,3),  yielding, 
o  _ 

G^  (e^,AT,a£)  —  ^Ic  (2.10) 

If  aT  Is  given,  it  follows  that  the  critical  (e^x^cr  formation  of  the 
first  transverse  crack  Is  determined  from  (2,10)  by  setting  a^  =  a.naici  yleld- 
ing 


^®x)cr 


/Gj({/2t  -  /Ct  •  A'f 

7^ 


(2.11) 


The  Propagative  Nature  of  Transverse  Crack,  in  order  to  see  the  propaga¬ 
tive  behavior  of  a  flaw  from  Its  Initial  size  23^  to  a  full  transverse  crack, 

o 

It  Is  useful  to  examine  the  behavior  of  G  In  (2.9),  which  Is  a  function  of  a. 

o 

Figure  2.10  shows  the  mode-I  energy  release  rats,  G  ,  plotted  against  a, 
typically  for  [O/OOla  laminates  under  uniaxial  tension. 

It  Is  seen  that  the  available  crack-tip  energy  release  rate  Is  a  non¬ 
linear  function  of  a.  The  general  behavior  Is  that  It  Increases  sharply  from 
Initial  value  of  aj  and  reaches  a  maximum  at  about  3/4  the  hal^  thickness  of 
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the  90‘*-plle8  (-  nt) .  After  this  point,  a  decrease  In  Is  effected  as  the 
crack  propagates  near  the  0/90  Interface.  Presumably,  the  crack  is  then 
arrested  at  the  Interface. 

If  the  Initial  size  aj  Is  sinaller  than  3/4  of  (nt) ,  onset  of  propagation 

o 

Is  determined  by  equating  G  (a^)  to  according  to  (2.10).  When  the  flaw 

o 

size  starts  to  Increase  from  a^,  the  available  G  becomes  larger  than  G^^.. 
Consequently,  the  propagation  la  unstable,  according  to  the  stability  criter¬ 
ion  of  (1.4). 

The  excess  energy  (the  shaded  area  In  Figure  2.10)  released  during  the 

unstable  propagation  caay  result  In  local  damages  at  the  0/90  Interface  near 

the  crack  termini.  This  subject  will  be  further  discussed  In  a  later  chapter. 

But,  It  Is  noted  that  the  amount  of  the  excess  energy  Increases  with  the 

number  of  the  OO^-plles,  n.  It  may  thus  be  Inferred  that  the  level  of  damage 

following  the  transverse  crack  arrestment  at  the  0/90  Interface  Increases  with 

the  number  of  the  OO’-plles  grouped  together. 

Simulation  of  Multiple  Cracks.  In  the  foregoing,  the  conditions  for 

determining  the  critical  onset  load  for  the  first  transverse  crack  have  been 

established.  After  the  forinatlon  of  the  first  crack,  all  future  cracks  will 

form  froin  the  remaining  flaws.  But,  the  sl'ies  of  the  remaining  Flaws  are 

smaller  than  the  first,  and  some  of  the  flaws  'nay  now  he  under  the  Influence 

o 

of  the  shear-lag  effect  from  the  First  crack.  Now,  let  G  be  the  available 

T 

energy  release  rate  at  the  tlu  of  a  flaw  If  It  Is  outslle  the  Influence  of 

shenr-lag,  while  G^  he  the  energy  releT<?e  rnte  If  It  Is  Inslle  the  shenr-lag 

o 

Then,  the  value  of  Is  smaller  than  G  and  Gj  Is  exi)resslhle  la  the 


zone. 
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where  R(s)  is  a  reduction  factor  representing  the  effect  of  shear-lag.  Depen¬ 
ding  on  the  distance  s  between  the  flaw  and  the  first  crack,  R(s)  has  a  value 
In  the  Interval  (0,  1),  see  Figure  2.11. 

The  curve  shown  In  Figure  2.11  can  be  numerically  generated  by  the  same 
finite  element  procedure  by  Including  the  presence  of  a  transverse  crack  In 
the  calculation  of  G  for  the  considered  flaw;  see  inset  of  Figure  2.11, 

Hence,  for  all  the  remaining  flaws  In  {a^},  minus  a,iiax»  flaw  that 
gives  the  largest  will  be  the  next  to  propagate  into  a  transverse  crack. 
This  Is  accomplished  most  routinely  by  a  computer  search  [12].  Then,  the 
corresponding  applied  laminate  tension,  e^,  for  this  flaw  to  propagate  is 
determined  again  from  the  fracture  condition, 

“  ^IG  (2.11) 

given  the  values  of  AT  and  a  associated  with  the  flaw. 

Note  that  determined  from  (2.11)  is  larger  than  (e'x)cr  determined  for 
the  first  crack,  because  a  Is  smaller  than  a^^^  and  the  flaw  may  be  <inder  the 
Influence  of  the  shear-lag  of  the  first  crack. 

As  for  subsequent  cracks,  they  can  be  similarly  determined  from  the 
remaining  flaws.  The  corresponding  laminate  strains  are  found  through  solving 
equations  similar  to  (2.11).  However,  a  flaw  located  between  two  adjacent 
transverse  cracks,  will  be  subjected  to  a  shear-lag  effect  from  each  crack. 
In  such  cases,  the  available  energy  release  rate  at  the  flaw  tip  Is  expressed 
as 

Gi(ex,  AT,  a)  =  R(s,)  •  G°  (I^,  AT,  a)  •  R(sr)  (2.14) 

where  S|^  and  Sj.  are  distances  from  the  flaw  to  the  left  crack  and  to  the  right 
crack,  respectively. 
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Successive  searches  for  the  most  eaergetlcally  possible  flaws  to  become 
transverse  cracks  provide  a  load-sequence  for  the  multiple  cracks  formation 
process.  The  searching  procedure  mlmlcks  the  cracking  process  as  It  woulo 
occur  naturally.  Since  each  simulation  begins  by  generating  a  set  of  M  random 
numbers  used  to  define  the  flaws  In  the  90®-plles,  the  random  number  set  may 
be  considered  In  a  sense  as  a  test  specimen.  Generations  of  several  random 
number  sets  can  be  used  to  simulate  replicate  test  specimens.  In  this  manner, 
the  scattering  property  of  transverse  cracks  In  a  sample  of  replicates  can 
also  be  simulated. 

2.6  Illustrative  Examples 

With  the  modelling  procedure  outlined  In  the  preceding  section,  the 
method  can  now  be  Illustrated  by  several  experimental  examples.  ConsHer 
first  the  T300/934  laminate  family  [0/90,^/0l ,  n  =  1,  2,  3  and  4  under  uniaxial 
tension.  The  basic  ply  properties  of  the  unidirectional  system  have  been 
previously  given  In  Table  TT-1,  and  the  growth  data  for  multiple  cracks  under 
Increasing  laminate  tension  were  shown  earlier  in  Figure  2.6.  From  this  crack 
growth  data,  the  experimental  onset  load  for  transverse  cracking  for  each 
value  of  n  can  be  determined.  These  are  listed  In  Table  II-2. 

Detarminatton  of  Effective  Flaw  Distributions.  The  first  step  In  the 
modelling  Is  to  determine  the  effective  flaws  In  the  basic  bO°-ply  as  a  ply 
property.  To  do  so,  a  correlation  la  performed  between  the  predicted  and  the 
experimental  crack-density  versus  applied  load  relationships.  ^rom  Figure 
2.6,  It  is  noted  that  In  the  case  of  n  =  l  there  was  no  noticeable  transverse 
cracks  formed  prior  to  laminate  final  rupture.  but  For  n  =  2,  3  and  4, 
sufficient  number  of  transverse  cracks  are  recorded.  For  the  llbistratlve 
examples  to  be  discussed  In  this  section,  the  effective  flaw  distribution  la 
the  90°-pl/  will  he  determined  using  the  data  curve  for  u  =  2. 
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assumed  here  that  two  (2)  90‘’“plles  grouped  together  shall  be  a  basic  unit; 
the  flaws  In  Its  thickness  are  denoted  by  the  set  {(32)1}.  Then,  according  to 
the  volume-law  proposed  In  Equation  (2.8),  the  flaw  sizes  for  n  =  3  and  n  =  A 
are  represented  by 


2/a- 

(an)i  =  (32)1  (n/2)  2 


3,  4 


(2.1‘>) 


where  a2  Is  the  shape  parameter  of  the  distribution  function  for  {32). 

Now,  assume  there  are  M  flaws  In  each  unit  length  of  the  basic  90°-ply; 
and  there  are  M  locations  (xi)  In  a  unit  length.  The  set  of  flaw  locations, 
{xi).  Is  obtained  from  a  suitable  random  number  generator  that  yields  M  values 
In  the  Interval  (0,  1).  These  same  values  are  also  used  In  conjunction  with 
the  assumed  flaw  distribution  function  of  the  Welbull  form.  Equation  (2.1),  in 
the  following  way;  by  letting  F  equal  to  each  of  the  discrete  values  of  (xi) 
and  by  guessing  the  parameters  02  and  S2  (the  subscript  2  denotes  two  90°- 
plles  grouped  together),  a  discrete  set  {(32)1} »  1  *  1,  M,  Is  obtained  from 
solving  (2.1). 

With  {xj}  and  {(32)1}  defined,  simulation  of  the  transverse  crack  process 
may  begin.  Detailed  calculations  In  the  simulation  will  be  presented  later  In 
this  section;  but  It  Is  sufficient  to  state  here  that  correlations  between  the 
simulated  results  and  the  experimental  data  finalize  the  values  of  a?  and  p? 
and  the  value  of  M.  For  the  data  shown  for  n  =  2  In  Figure  2.6,  this  proce¬ 
dure  yielded,  by  best  fit,  for  the  following  values: 


N  =  32  crack/cm 
12  =  4;  p,2  =  0,96t 


(2.16) 


Now,  let  the  largest  flaw  la  {(an)!!  be  In  the  top  3%  rank,  or  the  97- 
percentlle  In  the  flaw  size  cumulative  function  (2.1);  It  follows  that 
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n 


I 


(®2^max  ^  1.31t  (2.17) 

Using  the  volume-law  of  (2.15),  the  maximum  flaw  size  for  the  cases  of 
3  and  n  =  4  are  represented  respectively  by: 


(®3^max  ^  1.61t 
(®4^max  1.35t 


(2.18) 


Determination  of  Onset  of  Transverse  Cracking.  By  using  the  ply  stiff¬ 
ness  properties  of  the  T300/934  graphlte/epoxy  system  listed  In  T^ble  TT-1, 
the  strain  energy  release  rate  coefficients  Cg  and  Op  that  appear  In  Equation 
(2.9)  are  generated  numerically  for  each  value  of  n.  These  coefficients  are 
shown,  respectively,  in  Figures  2.12  and  2.13  as  functions  of  the  half  flaw 
size  a.  To  predict  the  critical  tensile  load  for  onset,  Equation  (2.11)  Is 
used  In  which  the  coefficients  Cg  and  Op  take  their  values  at  a  =*  '’max*  These 
values  are  listed  In  Table  11-3.  Note  that  for  the  case  of  n  =  1,  ’max  ould 
not  be  determined;  so  the  maximum  available  energy  release  rate  used 
for  prediction. 

finally,  by  assigning  2t  =  132  x  10”^m,  aT  =  IZS^C  and  r,i^  =  228  J/m^ 
(Table  TI-l)  solving  Equation  (2.11)  yields  the  critical  strain 
(ex)cr  each  case  of  n.  The  predicted  results  and  those  found  experiment¬ 
ally  (see  Table  TI-2)  are  compared  below: 


n 

^-x^cr*  Predicted 

(®x^cr»  Experiment 

.n 

1.05’i 

(^  1.0%) 

2 

0.68% 

0.66'’' 

3 

0.53% 

0.61% 

4 

0.52% 

0.52% 
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For  all  practical  purposes,  the  prediction  Is  satisfactory  both  In  njagnl- 
tude  and  In  trend. 

Simulation  of  Multiple  Cracks.  Since  the  results  shown  In  Figure  2.6  are 
from  the  Initial  portion  of  the  transverse  cracking  process.  It  Is  Insuf¬ 
ficient  for  a  complete  simulation.  Instead,  the  eKperlmental  data  reported 
recently  in  Reference  [38]  will  be  used  next.  In  that  experiment,  a  series  of 
T300/934  graphite-epoxy  laminates  In  the  form  of  [02/902n/^2l »  n  =  1,  2  and  4, 
were  tested  under  uniaxial  tension.  This  family  of  laminates  Is  Identical  In 
stacking  sequence  as  the  t0/90ij/0l  family  considered  In  the  earlier  examples, 
but  the  thicknesses  of  the  layers  of  like  fiber  orientation  are  now  doubled. 
As  It  will  be  shown  shortly,  the  doubling  of  the  layer  thickness  permits  early 
onset  of  transverse  cracks  In  the  90“-plles.  Consequently,  complete  develop¬ 
ment  of  transverse  cracks  Is  obtained  before  the  other  complicating  failure 
modes  occur. 

Figure  2.14  shows  the  experimental  crack  density  versus  the  applied  lami¬ 
nate  stress  Ox  curves  for  each  n  In  the  (02/^'^2n^^2l  family.  For  each  n, 
there  were  four  replicate  specimens  tasted;  and  a  finite-width  scatter  band  Is 
obtained  and  displayed.  Note  that  for  each  n,  a  transverse  crack  density 
saturation  limit  Is  reached  before  final  laminate  failure.  The  random  search 
procedure  outlined  earlier  will  now  be  used  to  simulate  this  result. 

The  T300/934  Graphite-epoxy  ply  constants  listed  In  Table  lT-1  are  used 

In  all  computations.  In  fact,  the  energy  release  rate  curves  for  Cg  and  Cp 

shown  In  Figures  2.12  and  2.13  are  also  used  In  the  present  [02/907ri/02] 

cases.  Since  layers  of  like  fiber  orientation  are  now  doubled,  the  curves  in 

o 

Figures  2.12  and  2,13  can  be  used  to  compute  G  In  equation  (2.9)  by  replacing 
2t  by  4t.  Then,  all  other  equations  remain  unchanged,  except  the  slv-e  distri¬ 
bution  function  for  the  effective  flaws. 
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As  the  effective  flaws  In  the  basic  piles  (two  90®-plles  grouped  toge¬ 
ther)  were  already  determined,  see  (2.16),  It  Is  necessary  only  to  follow  the 
volume-law  to  define  the  flaw  sizes  In  90**-layer8  of  more  than  two  piles. 

Note  that  In  this  laminate  series,  with  n  -  1,  2,  and  .4,  the  numbers  of 
the  grouped  90“-pllea  are,  respectively,  2,  4,  and  8.  The  effective  flaw  size 
for  the  cases  of  n  ••  1,  2  and  4  should  now  be  determined  by, 

2/  a« 

^®2n)l  *  (®2)l  ^  n  “  1,  2  and  4  (2.19) 

It  will  be  assumed  that  for  each  case  of  n,  there  are  M  »  32  flaws  per 
1  cm  as  before.  The  flaw  locations  are  again  represented  by  a  set  of  M  random 
values  In  the  Interval  (0,  1).  To  Include  the  shear-lag  effect,  the  energy 
reduction  factor  ^(s)  must  be  generated  Independently.  The  R(s)  factor  for 
I02/902n/02l  represented  by  a  single  function  In  terms  of  s  and  n: 

^(s,  n)  =  1  -  1.0056  exp  (-0.3838s/n  t]  (2.20) 

Equation  (2.20)  Is  graphically  displayed  In  Figure  2.15. 

Now,  all  the  elements  for  the  simulation  are  now  given:  the  value  of  M, 
the  flaw  size  distribution  {aj}  for  any  value  of  n,  the  flaw  location  distri¬ 
bution  {xj},  the  energy  release  rate  coefficients  Cg,  and  Op  and  the  energy 
reduction  factor  R(8,  n) .  By  following  the  outlined  search  procedure,  90®- 
layer  crack  development  Is  simulated  as  a  function  of  ascending  laminate 
stress  ax»  or  strain  e^.  The  computer  routine  utilized  Incorporates  a  Monte- 
Carlo  type  search  procedure,  the  detail  of  which  has  been  described  In  [I2l. 

For  each  [02/902n/02l  laminate  type,  four  or  five  replicate  simulations 
are  performed  (each  simulation  represens  n  specimen).  Because  of  the  discrete 
nature  of  the  random  numbers,  specimens  so  represented  are  slightly  different 
from  one  another.  Consequently,  repeated  simulations  also  give  rise  to  a 
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scatter  band  In  the  results. 

Figure  2.16  shows  the  simulated  transverse  crack,  density  versus  the 
applied  (tensile)  laminate  stress  ox  fot  four  specimens  of  n  =  1.  It  Is  seen 
that  the  simulated  results  form  a  narrow  scatter  band  and  agree  closely  with 
the  band  obtained  by  actual  experiment. 

Similar  comparisons  for  the  laminates  of  n  =»  2  and  4  are  displayed  In 
Figure  2.17  and  Figure  2.18,  respectively.  Here,  again,  the  simulated  results 
agree  reasonably  well  with  the  experiment  both  In  trend  and  In  scatter.  Note 
that  the  scatter  band  In  the  simulated  results  stems  from  the  random  numbers 
scheme  used,  whereas  the  scatter  band  In  the  experimental  data  Is  due  to 
physical  factors  such  as  the  random  property  of  the  basic  composite  material 
system  and  the  randomness  in  the  laminate  curing  procedure,  specimen  handling 
practice,  the  test  method,  etc.  Tn  the  Illustrations  shown  in  Figures  2.16  to 
2.18,  no  special  effort  was  made  to  correlate  In  any  way  the  scatter  between 
the  simulated  and  the  experimental  results. 

2 . 7  Effects  of  Thermal  kesldual  Stresses 

In  the  example  problems  considered  so  far,  the  effects  of  thermal  resid¬ 
ual  stresses  liave  been  Included  implicitly  In  the  analysis.  It  Is  Interesting 
to  examine  these  effects  explicitly  in  some  detail.  Figure  2.19  shows  the  In- 
sltii  tensile  stress  Ox  li  the  90®-plles  at  the  onset  of  the  first  transverse 
crack  for  each  n  value  In  the  [0/90i,/')]  family.  The  calculated  stress 
contains  the  contributions  from  both  the  thermal  cooling  of  M  and  the  mechan¬ 
ical  loading,  ex* 

It  Is  seen  that  the  total  In-sltu  ox  Is  not  constant  with  respect  to  n, 
but  decreases  sharply  with  Increasing  n.  Although  the  In-sltu  ox  1'’  directly 
responsible  for  the  formation  of  transverse  cracking.  It  Is  not  a  useful 
strength  measure  for  the  hO^-ply.  Thermal  cool  down  following  fabrication 
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curtag  Induces  about  25%  to  40%  of  the  total  In-sltu  Ox»  depending  on  the 
value  of  n. 

The  thermal  stress  contribution  to  the  total  strain  energy  release  rate 

o 

is  associated  with  the  coefficients  Cf  In  Equation  (2.9).  Note  that  0^  Is 

composed  of  three  parts:  Ggi  G-p  and  GeT  representing,  respectively,  the 

energy  due  directly  to  the  applied  tension  temperature  drop  AT  and  the 

coupled  energy  Involving  both  and  aT.  Figure  2.20  displays  the  values  of 

the  parts  Gg,  Gj  and  Gg'f  for  tlie  considered  laminate  family  when  the  first 

o 

crack  la  formed.  Note  that,  at  onset  of  the  crack,  G^  Is  equal  to  Giq  and  Gj(] 
is  Independent  of  n. 

Since  the  sum  of  G'l*  and  Gg-p  Is  attributable  to  the  thermal  effects.  Its 
share  In  the  total  energy  release  rate  can  be  as  much  as  40%  to  65%,  depending 
on  the  value  of  n. 

Clearly,  the  theraial  effects  are  Important  and  must  be  Included  In  the 
cracking  analysis.  Furthermore,  the  thermal  stress  effect  and  the  applied 
stress  effect  on  the  total  strain  energy  required  for  cracking  are  not  simply 
additive.  This  Is  fundamentally  different  from  the  failure  criterion  based  on 
strength.  In  the  latter,  thermal  and  mechanical  stresses  are  linearly  addi¬ 
tive. 


2.8  Summary 

In  this  chapter,  an  analysis  method  was  presented  for  the  formation 
process  of  Intraply  matrix  cracks  la  the  form  of  90°-ply  transverse  cracking. 

The  underlining  assumption  of  the  method  Is  that  successive  cracks  In  the 
90®-ply  under  transverse  tension  are  due  to  the  propagation  and  the  arrestment 
of  material  flaws  which  distribute  randomly  In  the  90®-ply.  The  analysis 
entails  two  essential  elements,  the  first  being  the  proper  representation  of 
the  material  flaws  as  an  Inherent  ply  property,  and  the  second  being  the 
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establishment  of  a  logical  criterion  determining  the  exact  loading  condition 
at  which  a  particular  flaw  becomes  a  transverse  crack. 

The  representation  of  material  flaws  Is  based  on  a  set  of  assumed  effect¬ 
ive  flaws  for  the  basic  material  ply.  As  such,  the  effective  flaws  do  not 
have  real  physical  Identity  but  their  effect  results  In  the  macroscopic 
characteristics  of  the  multiple  cracking  observed  In  experiment. 

Once  the  effective  flaw  distribution  for  a  90®-ply  having  a  reference 
thickness  Is  known,  the  effective  flaws  In  several  90®-plle8  grouped  together 
can  be  determined  by  means  of  a  statistical  volume-law.  The  law  suggested  by 
Equation  (2.8),  while  plausible.  Is  essentially  empirical.  Further  research 
Is  needed  In  order  to  formulate  a  volume— law  on  more  rigorous  grounds. 

Ka  for  the  formation  of  transverse  cracks,  the  classical  fracture  mechan¬ 
ics  criterion  Is  employed.  It  Is  assumed  that  each  of  the  effective  flaws 
behaves  like  a  small  crack.  Onset  of  Its  propagation  defines  the  transition 
from  a  flaw  to  a  transverse  crack;  and  Its  arrestment  at  the  0/90  Interface 
completes  the  formation. 

The  calculation  of  the  strain  energy  release  rate  associated  with  a  flaw/ 
crack  propagating  In  the  thickness  of  the  hO^-plles  Is  accomplished  by  a 
finite  element  procedure  based  on  ply  elasticity.  The  procedure  Is  a  numeri¬ 
cal  alternative  for  problems  having  complex  geometrical  constraints,  along 
with  other  peculiarities  such  as  the  Inclusion  of  thermal  effect,  shear-lag 
effect,  etc. 

Because  the  effective  flaws  In  the  90®-plles  are  random  In  sl^.e  and 
location,  successive  flaw/crack  transitions  under  the  applied  load  constitute 
a  stochastic  process.  ^  search  procedure  Is  used  for  the  simulation  of 
multiple  cracks  as  a  function  of  continued  loading. 


• h •  . 
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In  addition  to  the  logical  development  of  the  analytical  model,  the 
Influence  of  the  lamination  parameters  on  the  model  are  also  addressed.  These 
parameters  are  found  to  play  an  Important  role  In  the  crack  formation  process¬ 
es.  For  example,  the  effects  of  lamination  stacking  sequence,  the  actual  ply 
thickness,  etc.,  are  Implicitly  contained  In  the  predictive  model.  By  varying 
these  parameters  both  In  the  model  simulation  and  In  the  experiments,  their 
Influence  can  be  quantified  and  the  reliability  as  well  as  the  generality  of 
the  analytical  model  can  be  evaluated. 

At  the  conclusion  of  this  chapter,  it  Is  noted  that  the  effective  flaw 
concept  results  In  physically  correct  behavior  although  the  actual  representa¬ 
tion  of  the  flaws  remains  a  subject  for  further  refinement.  This  also 
Includes  the  empirically  derived  volume-law  for  flaws  in  grouped  90“-plles. 
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TABLE  II-l 

Material  Constants  for  T300/934  Unidirectional  Ply* 


Property 


SI  Unit 


English  Unit 


Ell 

144.8  Gpa 

21x10^  psl 

Err*  Ezz 

11.7  Gpa 

1.7x10^  psl 

VLTt  VLz 

0.3 

0.3 

VTz 

0.54 

0.54 

GlT»  ^Lz 

6.5  Gpa 

0.94x10^  psl 

^Tz 

3.5  Gpa 

0.5x10^  psl 

OT*  Oz 

?8.8x10“6/“C 

16.10~6/“F 

aL 

0.36.10-^/“C 

0.2x10"6/“F 

OL,  long,  tensile  strength 

1720  mpa 

250  ksl 

ox»  trans.  tensile  strength 

44  mpa 

6.5  ksl 

2t,  nominal  ply  thickness 

0.132  mm 

0.0052" 

at,  temperature  drop 

125‘’C 

225'’F 

Gic  (0/0  Interface)** 

158  J/m2 

0.9  In-lh/ln^ 

Gic  (90/90  Interface)** 

228  J/'n2 

1.3  In-lb/ln^ 

G(I,TI)c»  =  9.5** 

288.75  .J/n2 

1.65  In-lb/ln^ 
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TABLE  II-2 

Onset  of  Transverse  Cracking  In  [0/90f|/0l  Laminate  Faintly* 


Lamina te 


Experimental 
Axial  Stiffness 


Experimental 
Onset  Stress 


Experimental 
Onset  Strain 


[0/90/0] 

95.8  GPa 

>950  MPa 

>1.00% 

[0/902/0] 

78.2  GPa 

520  MPa 

0.66% 

[O/9O3/O] 

64.9  GPa 

400  MPa 

0.61% 

[O/9O4/O] 

55.8  GPa 

290  MPa 

0.52% 

*Dat,a  taken  from  Reference  i32l 
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TABLE  II-3 

Maximum  Flaw  Size  and  Energy  Release  Rate  for  [0/ 


n  -e.  cT,  l/.n3 


♦Maximum  value  Is  used. 
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Figure  2.1  The  Reinforcement  Structure  of  Laminates  Made  of  Unidirectional  Plies 
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d/t  16  12  8  4 


Figure  2.3  Interlaminar  (in  kpa/ye)  on  Laminate  Mid-Plane 
as  a  Function  of  Distance  From  Free  Edge  [0/90] g. 


Figure  2.6  Experimental  Crack-Density  versus  Laminate  Tension  Relations  for  the 
[0/90jj/0]  Laminate  Family. 
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Figure  2.7  An  X-Radiograph  and  a  Thickness  View  Photo-Micrograph 
of  Transverse  Cracks  in  a  [+25/902] g  Laminate. 
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transv.  crack 


Excess  Energy  Release  Rate  Given  Off  by  a  Transverse  Crack  PropagatinB  In 
Thick  90®-layer. 


Figure  2.11  Energy  Release  Rate  Retention  Factor  R(S)  for  the 
[0/90^/01  Laminate  Family. 
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10*  J/m* 


0  1  2  3  4  a/t 


Figure  2.12  Energy  Release  Rate  Coefficient  Cg  for  the 
(0/90jj/0]  Laminate  Family. 
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ck8/cm 


Experimental  Crack-Density  versus  Lamination  Relations  for  the  [02/90^] 
Laminate  Family. 


.15  Energy  Release  Rate  Retention  Factor  R(S)  for  the  [02/90^13  Laminate  Family. 
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Laminate  Tension  Relations  [O2/9O] 


SIMULATION 
4  SPECIMENS 
C  02/904  3s 
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Figure  2,18  Comparison  Between  the  Simulated  and  the  Experimental  Crack-Density 
versus  Laminate  Tension  Relations  [02/904]s. 
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CHAPTER  lit.  ANALYSIS  OF  INTERPLY  CRACKING 

3,1  Physical  Charactertsttcs  of  Interply  Cracking 

In  Chapter  II,  the  oiechanlcs  of  Intraply  cracking  In  the  form  of  trans¬ 
verse  cracks  was  Investigated.  Intraply  cracking  Is  labeled  as  the  first 
basic  form  of  matrix  crack,  which  occurs  primarily  because  of  the  In-ply 
stresses,  Oy*  Oy,  and  T^y.  The  second  basic  form  of  matrix  crack  Is  Interply 
crack,  or  delamlnatlon.  Oelamlnatlon  Is  a  plane  crack  which  Initiates  and 
propagates  In  a  plane  parallel  to  or  Inside  of  an  Interface  between  two 
adjacent  piles.  The  primary  cause  Is  the  Interlaminar  stresses  t^z 

Ty2  that  exist  In  the  delaminating  plane. 

In  practice,  most  structural  laminates  are  designed  to  withstand  In-plane 
loads.  Interlaminar  stresses  occur  near  locations  of  geometric  discontinuity, 
such  as  laminate  free  edge,  bolt-hole,  surface  cut,  Inclusion,  Internal 
defect,  or  some  Intraply  cracks  that  already  occurred,  etc.  Hence,  Interlami¬ 
nar  stresses  ar?  localized  or  concentrated  stresses. 

The  severity  of  the  various  localized  Interlaminar  stresses  depends  on 
the  nature  of  t)ie  particular  geometrical  discontinuity.  Similarly,  the 
Initiation  and  propagation  characteristics  of  Interply  cracking  also  varv  with 
the  same  geometrical  factor. 

In  this  section,  an  effort  Is  inale  to  delineate  the  physical  behaviors  of 
some  commonly  observed  de lamina tlons  In  engl  leerlng  terms,  and  to  set  a 
logical  procedure  for  an  analytical  treatment. 

Free  Edge  Oelamlnattons.  Consider  first  the  Free  edge  delamlnatlon 
problem.  This  problem  has  attracted  reseircii  literests  ever  sloce  the  advent 
of  composite  laminates.  As  far  back  as  In  the  late  HGTs,  labor  itory  coupon 
tests  already  established  that  the  tanslle  strength  of  a  Timlmt’  coupon 
depends  on  whether  or  not  free  edge  delamlnatlon  occurs  prior  to  Final  coupon 
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rupture.  In  this  regard,  It  was  found  that  the  laminate's  ply  stacking 
sequence  greatly  Influences  the  behavior  of  the  observed  free  edge  delaralna- 
tlon  [39].  Bjeletlch,  et.  al  [29]  showed,  for  example,  the  uniaxial  tensile 
strength  of  the  graphite-epoxy  [+45/0/90] g  coupon  was  much  lower  than  the 
[0/90/+45]g  coupon.  The  reason  Is  that  the  former  can  develop  extensive  free 
edge  delamlnatlon  prior  to  final  failure,  while  the  latter  cannot.  Free  edge 
stress  analysis  based  on  ply  elasticity  reveals  that  large  tensile  stress 
la  concentrated  along  the  free  edge  of  the  [+45/0/90] g  coupon,  while  a 
compressive  stress  Is  concentrated  along  the  edge  of  the  [0/90/+45]3 

coupon.  Hence,  the  former  Is  prone  to  delamlnatlon;  the  latter  Is  not. 

Actually,  stacking  sequence  Is  not  the  only  factor  affecting  the  behavior 
of  free  edge  delamlnatlon.  Other  factors  liave  been  identified.  For  Instance, 
Rodlnl,  et.  al  [40]  were  the  first  to  experimentally  document  the  dependence 
of  delamlnatlon  growth  on  ply  thickness.  They  showed  that  the  critical 
laminate  tensile  strain  e^  at  the  initiation  of  edge  delamlnatlon  In  a 
[_+45^j/0n/90ji]3,  n  =  l,  2,  3,  laminate  family  varied  approximately  Inversely 
with  the  value  of  »'n.  A  more  systematic  study  of  ply  thickness  effect  on 
matrix  cracking  In  general  was  reported  recently  by  Crossman,  et.  al  [41]. 

Experimental  documentation  of  free  edge  delamlnatlon  has  traditionally 
been  conducted  using  both  macroscopic  and  microscopic  damage  Inspection 
techniques.  Optic  opaque  penetrant  enhanced  x-radlogr iph/  I-  nong  the 
techniques  used  at  the  macroscale,  while  microphotography  In  conjunction  with 
scanning  electronic  microscopy  Is  a  method  for  detecting  material  damage  at 
the  microscale. 

Figure  3.1a  shows  an  x-radlograph  In  plane  view  of  free  edge  delamlnatlon 
In  a  [^45/9/90]g  graphlte-epoxy  coupon  subjected  to  uniaxial  tension.  Tn  this 
particular  laminate  coupon,  a  uniform  delamlnatlon  was  Initiated  along  both 
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edges  of  the  laminate  at  about  90%  the  laminate  ultimate  strength.  The  growth 
behavior  was  essentially  stable.  Near  total  de lamination,  other  matrix  crack¬ 
ing  modes  developed  and  caused  final  failure  of  the  test  coupon. 

Microscopic  examination  of  the  same  delaminated  coupon  gives  a  more 
detailed  view.  Figure  3.1b  Is  an  edge  view  of  the  delamlnatlon  magnified  at 
80  times.  This  picture  was  taken  after  the  onset  of  edge  delamlnatlon,  but 
before  the  laminate  failure.  It  Is  seen  that  the  crack  la  not  confined  In  any 
one  given  ply-to-ply  Interface.  Rather,  It  wanders  In  the  thickness  of  the 
OO^-plles  with  a  xlg-zagged  path.  The  zig-zagged  crack  surface  may  be 
explained  by  the  fact  that  the  OO^-plles  near  the  free  edge  region  <xe  subi act¬ 
ed  both  to  a  tensile  normal  to  the  laminate  plane  and  to  a  tensile  ox 
the  loading  direction.  The  tensile  Ox  tends  to  cause  transverse  cracks  In  the 
OO^-ply,  while  the  tensile  oz  causes  delamlnatlon. 

Ka  the  applied  tensile  load  Increases,  the  delamlnatlon  propagates 
further  toward  the  center  of  the  laminate.  The  growth  Is  sometimes  concurrent 
with  multiple  transverse  cracks  In  the  90°-ply.  i^'lgure  3.3  shows  an  80x 
thlckness-v lew  of  the  same  [_^45/0/90] 3  coupon  tsken  Just  before  Its  final 
failure.  Note  that  In  addition  to  delamlnatlon  In  the  hO^-plles,  there  are 
also  delamlnatlons  In  the  +45  Interfaces  and  the  -45/90  Interfaces.  Moreover, 
Intraply  cracks  In  the  +45  piles  are  also  visible.  Clearly,  at  this  stage  of 
loading,  laminate  matrix  cracking  Involves  several  modes  of  failure  that  occur 
at  several  Interfaces.  For  the  laminate  shown  here,  however,  onset  of  edge 
delamlnatlon  at  90%  of  the  laminate  ultimate  strength  appears  to  loe  the  first 
event  In  the  matrix  cracking  process. 

4s  mentioned,  growth  of  delamlnatlon  causes  isecondary  cracks  '0  other 
piles.  Figure  3.3  shows  a  sat  of  free  edge  delamlnatlon  phot)graphs  tikan 
from  a  ( +-25/90/ F25 j  tensile  coupon.  Mere,  a  free  elg.e  crack  ts  seen  to  have 
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formed  Initially  In  the  00“-ply,  but  then  branched  onto  the  25/90  Interface  as 
It  propagated  Inward.  The  branching  Is  probably  caused  by  the  secondary 
skewed  cracks  that  occur  ahead  of  the  delamlnatlon.  The  skewed  secondary 
cracks  In  the  90°-ply  thickness  Is  clearly  the  result  of  multl-axlal  tensile 
stresses  that  develop  near  the  delamlnatlon  front.  In  this  particular 
Instance,  the  thickness  of  the  90“-ply  Is  thin  and  the  threshold  condition  for 
a  pure  90°-ply  transverse  cracking  Is  not  reached  before  edge  delamlnatlon. 
However,  the  tensile  Ox  significant  enough  to  alter  the  Initial  delamlna¬ 
tlon  growth  direction. 

Finally,  Figure  3. A  displays  an  x-radlograph  of  free  edge  delamlnatlon 
for  a  [9O2/O2/+A52] s  laminate  coupon  subjected  to  uniaxial  compression  [42] . 
In  this  case,  a  tensile  stress  Is  developed  In  the  H^5°-plles  near  the  free 
edge.  This  caused  free  edge  delamlnatlon  In  the  A5/-A5  Interface.  Note  that 
the  delamlnatlon  contour  Is  not  as  uniform  along  the  Free  edge,  due  possibly 
to  the  constraint  from  the  grips  at  the  coupon’s  end-tabs.  Although  the  test 
specimen  was  very  short  length-wise  to  prevent  global  buckling,  local  ply 
buckling  still  occurred  subsequent  to  delamlnatlon.  The  latter  had  resulted 
In  transverse  cracks  In  the  A5°-plles,  as  can  be  seen  In  Figure  3.4. 

Other  Types  of  Oelamlnatlon.  The  foregoing  discussions  Focussed  on  the 
physical  nature  of  free  edge  delamlnatlons  as  observed  In  some  laboratory 
coupon  tests.  Generally  free  edge  delamlnatlon  Is  not  a  concern  la  applica¬ 
tions,  In  practical  problems,  delamlnatlons  usually  occur  near  laminate 
defects  where  the  stress  concentration  is  much  more  severe  than  the  Free  edge 
Induced  stress  concentration. 

figure  3.5  shows  an  x~radlograph  taken  From  a  graphl t^-?noxy 
( 45/O2/-A5/92/9OI  g  laminate  tensile  coupon  having,  a  throug.h  hole.  In  this 
test  coupon,  both  the  straight  edge  and  the  curv»el  ebie  aroimd  the  hoh»  caus.?! 
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interlaminar  stresses.  Indeed,  delamlnatlon  cracks  can  be  seen  to  form  along 
both  the  straight  edges  and  around  the  hole.  However,  the  contour  of  delaml¬ 
natlon  around  the  hole  Is  more  complicated  than  that  along  the  free  edge, 
because  the  stress  field  around  the  hole  Is  more  complex. 

Another  example  of  complicated  ply  delamlnatlon  Is  shown  In  Figure  1.6, 
where  a  T300/934  graphite-epoxy  tensile  coupon  having  [92/^*^3l3  lamination  was 
subjected  to  uniaxial  tension.  Purely  90“-ply  transverse  cracks  are  Induced 
long  before  any  delamlnatlon.  Upon  Increased  loading,  the  transverse  cracks 
acted  like  local  defects  which  cause  localized  Interlaminar  stresses.  This 
results  In  the  development  of  localized  delamlnatlons.  As  can  be  seen  from 
Figure  3.6,  these  delamlnatlons  grew  with  the  Increasing  load  and  their 
coalescence  precipitated  a  massive  0/90  Interface  separation.  This  last  event 
Is  believed  to  have  caused  the  final  specimen  failure. 

Finally,  Figure  3.7  Is  an  x-ray  picture  of  a  graphite-epoxy  laminate 
subjected  to  low  velocity  Impact.  The  localized  damage  caused  by  the  Impact¬ 
ing  object  Involved  primarily  matrix  cracks  In  the  various  plies  and  ply 
Interfaces.  Post-Impact  strength  Is  adversely  affected  by  the  damage  propaga¬ 
tion  upon  re-load Ing. 

These  laboratory  examples  are  cited  to  portray  the  complex  characters  of 
Interply  cracking.  Though  the  basic  cause  of  delamlnatlon  Is  attributable  to 
the  existence  of  Interlaminar  stresses,  the  exact  distribution  of  these 
stresses  cannot  always  be  determined.  Since  Interply  cracking  and  Intraply 
cracking  often  mutially  Interact,  the  Interacting  mechanisms  also  remain 
perplexing. 

In  the  following,  an  effort  Is  made  to  Idealize  and  simplify  the  lelaml- 
natlon  mechanisms,  yet  to  retain  all  the  es.sont  ini  Influencing  ractar,-.  Tn 
this  context,  a  descriptive  model  Is  formulated  for  the  Initiation  aal  growth 
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of  some  Important  delamlnation  problems,  as  seen  from  a  macroscopic  viewpoint. 

3.2  The  Energy  Model  for  Free  Edge  Delamlnation 

As  has  been  previously  stated,  the  physical  complexities  of  Interply 
cracking  are  far  from  being  fully  clarified,  especially  the  mechanisms  associ¬ 
ated  with  crack  Interactions.  But,  for  the  simpler  problem  of  free  edge 
delamlnation  which  propagates  without  crack  Interactions,  a  reasonably 
accurate  analysis  of  the  phenomena  at  the  macroscopic  level  Is  possible.  The 
analysis  approach  will  be  the  same  fracture  mechanics  concept  applied  to  the 
Intraply  cracking  problems.  Some  preliminary  considerations  for  the  model 
formulation  ace  given  In  the  following. 

Kinematics  of  Free  Edge  Delamlnation.  Consider  the  case  of  a  long  lami¬ 
nate  coupon  having  straight  edges.  Let  the  laminate  ply  stacking  sequence  and 
the  manner  of  loading  be  such  that  free  edge  delamlnation  occurs  exclusively 
at  a  known  ply  Interface.  For  simplicity,  assume  that  the  Initiation  and 
propagation  of  the  delamlnation  is  the  only  matrix  crack  event  during  the 
course  of  loading.  Under  such  circumstances,  the  delamlnation  growth  contour 
will  be  uniformly  advancing  along  the  edges;  the  growth  may  be  treated  as  one- 
dlmenslonal  and  self-similar.  This  Idealization  permits  the  representation  of 
the  delamlnation  as  a  line  crack  propagating  In  the  cross-sectional  plane  of 
the  laminate,  as  Illustrated  schematically  In  Figure  1.3. 

Effective  Interlaminar  Flaws.  For  definiteness,  let  the  laminate  be  of 
finite  width  2b  and  subjected  to  the  tinlaxlal  tensile  laminate  strain  e^. 
Then,  given  the  geometrical  lay-up  of  the  laminate,  all  free  edge  stresses 
Including  the  Interlaminar  ean  he  calculated  based  on  ply 

elasticity.  In  order  to  render  a  prediction  for  delamlnation  Initiation,  It 
Is  necessary  to  Introduce  a  certain  effective  Interface  flaw  near  the  free 
edge.  Upon  loading,  the  flaw  la  under  the  Influence  of  the  free  edge  stress 
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field.  Delamlaatlon  of  macroscopic  dimension  will  result  from  the  growth  of 
the  flaw  once  the  edge  stress  field  reaches  a  critical  state. 

Once  again,  the  size  of  the  effective  flaw  is  unknown;  some  heuristic 
arguments  about  micro-macro  transition  are  necessary  to  justify  Its  presence 
and  to  estimate  Its  size,  as  was  done  In  the  Intraply  cracking  problems. 

For  the  present  case,  let  It  be  assumed  that  an  "effective"  Interlaminar 
flaw  of  uniform  width  "Bq"  is  situated  In  some  ply  Interface  near  the  straight 
edge  of  the  laminate.  This  effective  flaw  Is  regarded  as  a  small  Interply 
crack  which  begins  propagation  when  the  stress  state  surrounding  the  flaw-tip 
reaches  the  critical  value.  This  Incident  then  defines  the  onset  of  de lamina¬ 
tion,  with  a  physical  crack  detectable  at  the  tnacroscale. 

Condition  for  Onset  of  Delamlnatlon.  The  transition  from  an  effective 
flaw  to  a  physical  crack  Is  governed  by  the  criterion  of  fracture  growth. 
Since  the  particular  interface  In  which  edge  delamlnatlon  occurs  is  assumed 
known,  then  the  flaw-tip  energy  release  rate  G(^x»  **o^  ^  given  effective 
flaw  of  size  a^  can  be  calculated.  The  finite  element  crack-closure  proce¬ 
dures  outlined  previously  Is  again  applicable.  The  calculated  G  can  be 
conveniently  expressed  as  a  function  of  ag,  given  the  stacking  sequence  of  the 
laminate,  the  ply  elastic  constants  and  the  applied  tension  e^. 

Figure  3.8  shows  a  typical  G-curve  plotted  against  the  flaw  size  ag.  It 
Is  seen  that  G  rises  sharply  from  zero  at  Sg  =  0  and  reaches  a  limiting  value 
Gmax  about  ag  >  The  fact  that  G  Initially  Increases  linearly  with  the 
flaw  size  ag  Is  a  behavior  similar  to  the  classical  crack  problems.  But  as  ig 
approaches  a^,  the  crack  driving  force  G  becomes  non-llnenr  and  reaches 
eventually  an  asymptote.  This  behavior  of  G  indicates  tiTat  the  delamlnatlon 
growth  Is  basically  a  stable  process  when  Its  size  exceeds  the  quantity  Sg. 


The  physical  meaning  of  Is  not  fully  understood  at  the  present  time. 
Its  numerical  value  for  most  laminates  Is  of  the  order  of  the  thickness  of  the 
material  piles  which  contains  the  delamlnatlon  crack.  For  example,  a^,  equals 
about  4t  (2  ply  thickness)  In  the  [+45/0/90] g  laminate;  It  equals  about  8t  In 
the  [+452/02/902]s  laminate  [32].  Hence,  the  value  of  similarly 

affected  by  the  same  ply  thickness  factor. 

The  calculated  G  function  can  be  expressed  again  in  explicit  terms  of  the 
applied  tension  e^t 

G(ex,  a)  =  Ce(a)  •  2t  •  e^^  (3.1) 

where  Cg(a)  Is  a  coefficient  dependent  only  on  the  delamlnatlon  size  (denoted 
as  a);  and  t  Is  a  chosen  length  parameter  (2t  =  ply  thickness).  Note  that  the 
functional  shape  of  Gg  Is  similar  to  G(a)  shown  In  Figure  3.8,  given  t  and 

Similarly,  If  the  laminate  Is  subjected  to  a  known  uniform  temperature 
drop,  at,  and  If  the  thermally  Induced  stresses  surrounding  the  flaw-tip  also 
effect  delamlnatlon  growth  In  similar  manner,  then  the  associated  flaw-tip 
energy  release  rate  due  to  aT  Is  expressible  as 

G(AT,  a)  =  Cx(a)  •  2t  •  AT^-  (3.2) 

where  Cq'(a)  Is  also  a  function  of  the  delamlnatlon  size  a  only.  The  function¬ 
al  behavior  07(3)  Is  similar  to  that  of  Cg,  see  (321. 

When  the  mechanical  load  and  the  temperature  load  both  effect  crack 
growth  at  the  same  time,  the  total  energy  release  rate  Is  expressible  is 

G(ex,  AT,  a)  =  [/Cg  •  e^  +  /^r  *  aTV?  (2t)  O.D 

Now,  if  the  Initial  size  of  the  effective  edge  flaw  a^  Is  V;nown,  t'u'ti  the 
condition  for  Its  propagation  (I.e.  onset  of  delamlnatlon)  is  governed  bv 
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G(®x»  ™  ^c» 


(3.4) 


In  order  to  execute  a  solution  for  (3.4),  one  question  must  now  be 
answered  and  two  quantities  defined.  The  question  Is  on  which  ply  Interface 
Is  the  delamlnaton  to  occur;  and  the  two  quantities  to  be  defined  are  the 
values  of  G^.  and  ag. 

The  Range  Of  Effective  Flaw  Size.  The  question  regarding  the  delamlna- 
tlon  site  and  the  value  of  will  be  examined  In  more  detail  latter  In  this 
section.  Now,  examine  first  the  possible  values  of  ag.  Note  from  Equation 
(3.4)  that  a  one-to-one  relationship  between  the  applied  tension  e^  and  the 
flaw  size  Sg  exists,  given  AT  and  Gg.  Figure  3.9  shows  a  schematic  plot  of 
Equation  (3.4)  In  the  e^  versus  flaw  size  plane.  It  Is  seen  that  If  a  range 
of  values  Is  given  for  the  edge  flaw  size  ag,  then  a  corresponding  range  for 
the  onset  tension  e^  can  he  determined.  In  particular.  If  for  all  probabili¬ 
ty,  the  effective  flaw  size  ag  Is  of  the  order  of  an,  or  larger,  then  the  rsnge 

of  the  onset  tension  e,^  Is  very  narrow.  Actually,  Equation  (3.4)  predicts  a 
constant  for  all  ag  >  a^,  because  G  then  reaches  Its  asymptote.  In  fact, 
this  constant  e-^  Is  the  minimum  possible  laminate  strain  predicted  for  which 
onset  of  delamlnatlon  could  occur. 

As  discussed  before,  the  numerical  value  for  a,n  In  many  practical  lami¬ 
nate  stacking  sequence  is  less  than  4t,  or  about  ?  ply  thicknesses.  It  Is 

entirely  possible  that  the  largest  affective  Interply  flaw  size  may  even  be 
greater  than  4t.  In  any  event,  a  lower  bond  for  the  onset  load  (e,^)gp.  can 
always  be  obtained  by  setting  ag  >  a^  or  Ug  ~  4t,  whichever  Is  smaller. 

Location  of  Delamlnatlon.  To  determine  the  exact  delamlnatlon  site  by  an 
analytical  means  Is  not  s tra Iglit-forward ,  because  given  a  multl-l  lyered 
laminate  coupon  there  are  always  several  ply  Interfaces  possible  for  delamtna- 
tlon. 


•  w  •  «  •  k  •  -  •  »  •  -  »  -  •  P  -.N  ^  •  .  •  .  •  .  N  . 
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Consider  for  example  the  T300/934  graphite-epoxy  (+25/90] g  laminate  under 
uniaxial  tension  examined  earlier  in  Chapter  II.  As  was  pointed  out  then, 
a  strong  tensile  Oz  is  concentrated  near  the  free  edge.  So,  free  edge  delaml- 
natlon  is  expected  before  OO^-ply  transverse  cracking.  The  interlaminar  oz 
distribution  along  the  mid-plane  (90/90  Interface)  was  shown  earlier  in  Figure 
2.4.  It  was  Inferred  that  it  is  this  stress  which  is  responsible  for  edge 
delamination. 

Actually,  Is  not  the  only  strong  interlaminar  stress  near  the  free 
edge;  nor  is  it  limited  to  act  only  on  the  laminate  mid-plane.  In  fact,  there 
is  a  high  shear  acting  both  on  the  2‘if~25  and  the  -25/90  interfaces  along 
with  a  tensile  Figure  3.10  shows  the  through-the-thlckness  distribution 

of  Txz  lesr  the  free  edge,  while  Figure  3.11  shows  the  through-the-thlckness 
distribution  for  oz«  is  seen  that  r^z  ^ce  singular  on  the  25/-25  and  the 

-25/90  interfaces.  Likewise,  a.  Is  singular  and  tensile  on  the  -25/90  inter¬ 
face. 

From  the  interlaminar  stress  distributions,  it  would  appear  that  the 
25/-25,  -25/90  and  90/90  Interfaces  are  all  possible  delamination  sites. 
Owing  to  the  lamination  symmetry,  delamlnatlon  la  the  90/90  interface  is 
predominantly  in  mode-I  (opening  mode);  delamlnatlon  la  the  ?5/-25  Interface 
is  in  mode-IIi  (anti-plane  shear)  and  delamlnatlon  la  the  -25/90  Interface  is 
mixed  with  mode-I  and  mode-Ill,  as  both  rxz  ^'id  act  on  that  Interface. 
From  such  a  stress  analysis,  however.  It  is  still  not  possible  to  determine 
for  sure  on  which  of  the  ply  Interfaces  delamlnatlon  will  initiate.  Clearly, 
a  suitable  quantity  Is  needed  to  explicitly  distinguish  the  delamiaatlon 
conditions  on  the  various  possible  Interfaces. 

To  this  end,  examine  the  computed  strain  energy  release  rate  avallahlie  la 
each  of  the  possible  interfaces.  dace  on  each  ply  interface*  an  effective 
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Interlaminar  flaw  of  size  ag  la  assumed,  the  available  total  G  at  the  flaw-tip 
in  each  of  the  interfaces  may  be  calculated  individually  in  terms  of  ag  and 
the  applied  tension  e^* 

Figure  3.12  shows  two  curves  for  the  total  energy  release  rate  coeffi¬ 
cient  Cg(a),  one  for  a  flaw  residing  in  the  mid-plane  (90/90  Interface)  and 
the  other  for  a  flaw  In  the  -25/90  Interface.  Note  that  Gg(a)  for  the  mid- 
plane  flaw  contains  exclusively  mode-I  energy,  while  CgCa)  for  the  flaw  In  the 
-25/90  Interface  (dotted  line)  contains  mode-I  and  mode-lll  energies.  The 
ratio  between  mode-IIl  and  mode-I  In  the  latter  Is  about  0.8. 

Similarly,  Figure  3.13  shows  the  energy  release  rate  coefficients  GT(a) 
for  delamlnatlon  In  the  same  two  Interfaces.  The  corresponding  Ce(a)  and 
G'r(a)  curves  for  delamlnatlon  in  the  25/-25  Interface  (predominantly  mode-III) 
are  omitted  because  of  their  relatively  small  magnitudes. 

From  Figures  3.12  and  3.13,  It  can  he  seen  that  propagation  of  a  flaw  In 
the  mid-plane  releases  more  energy  than  a  flaw  of  same  size  propagating  In  the 
-25/90  Interface.  However,  whether  or  not  a  mid-plane  delamlnatlon  will  occur 
before  a  -25/90  Interface  delamlnatlon  is  determined  by  the  specific  delamlna¬ 
tlon  resistances  Gjf]  for  mode-I  and  G(i^  llI)C  mixed-mode.  As  was  dis¬ 
cussed  In  Chapter  I,  different  cracking  modes  may  mean  different  microscopic 
failure  processes,  which  In  turn  give  rise  to  different  macroscop Ically 
measured  material  energy  release  rates. 

Dependence  of  Gg  on  Delamlnatlon  Modes.  In  several  recent  experiments  on 
graphite-epoxy  composites  [43  -  45l,  the  total  energy  release  rate  Gg  for 
mixed-mode  delamlnatlon  has  been  found  to  be  generally  greater  than  Gp-;  for 
mode-I  delamlnatlon.  In  fact,  In  one  case  the  total  Gg  for  ml<ed  mode  crack 
growth  was  found  to  Increase  with  the  ratio  G^j/Gj  [45] .  figure  3.14  shows 
the  measured  total  G~  for  mixed-mode  (I  and  TI)  delamlnatlon  plotted  as  a 


NADC-85118-60 


function  of  Gjj/Gj  ratio  (the  inateclal  used  was  the  T300/934  composite 
system).  It  shows  that  the  total  Is  always  greater  than  Gjq;  and  there 
seems  to  exist  a  limiting  value  at  large  Gjj/G-j.  This  limiting  value  may  be 
construed  as  the  value  of  Gx^q. 

The  reason  for  the  dependence  of  the  total  strain  energy  dissipated  in 
mixed-mode  crack  on  the  Gxi/Gx  ratio  may  be  explainable  at  the  microscale. 
Some  have  attributed  the  Increased  energy  to  crack-tip  matrix  yielding  under 
shearing  deformation  [44] ;  others  attempted  to  explain  It  from  the  crack 
surface  morphology  [43].  Whatever  the  exact  cause,  the  macroscopic  G^  Is 
multiple-valued  In  mixed-mode  crack  growth.  When  used  In  the  context  of 
Equation  (3.4),  It  becomes  necessary  to  know  the  Gxx/Gi  ratio  of  the  delamlna- 
tlon  growth,  so  as  to  select  a  proper  value  for  GQ-total.  Eor  this  reason,  a 
functional  curve  for  each  material  system,  such  as  shown  in  Figure  3.14,  must 
be  experimentally  determined. 

Alternative  mixed-mode  crack  growth  criteria  have  been  suggested  (see, 
e.g.  [46]),  requiring  material  resistances  such  as  Gxxc  ’’nd/or  Gxi;tC* 
many  questions  need  to  be  settled  concerning  both  the  existence  of  and  test 
methods  for  measuring  Gx[c  and  Gxjxf;  as  applied  to  de lamination  in  composite 
laminates.  For  the  time  being  the  criterion  In  Equation  (1.3)  will  be  used  in 
conjunction  with  a  variable  fracture  toughness  property,  such  as  displayed  In 
Figure  3.14  for  the  T300/934  system.  As  further  research  clarifies  the  'many 
uncertainties  surrounding  the  problem  of  nlxed-mode  matrix  cracking  In  compo¬ 
sites,  a  more  precise  and  direct  delamlnatlon  criterion  mny  be  developed. 

3.3  Illustra tlve  Example s  for  Fr ee  Ed ge  Oelamlna t Ion 

In  this  section,  two  example  problems  will  be  discussed.  All  are  related 


to  free  edge  delamlnatlon.  Hut,  the  loading  factor  and  lamination 
factors  v/lll  be  varied  to  see  how  rjood  or  bad  the  nredlctlve  model  Is 


geoiiie  try 
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[-25/90] 3  T300/934  Tenstoa  Coupoa.  Consider  the  [_+25/90]3  tensl  on  coupon 
made  of  T300/934  discussed  earlier.  A  prediction  can  now  be  made  as  to  which 
Interface  an  edge  delamlnatlon  will  first  Initiate.  If  the  fracture  criterion 
(3.4)  is  used  In  conjunction  with  the  total  G(,  curve  for  mixed-mode  delamlna¬ 
tlon  shown  In  Figure  3.14,  It  may  be  concluded  definitely  that  free  edge 
delamlnatlon  will  occur  In  the  mid-plane  of  the  laminate  coupon,  based  on  the 
calculated  energy  release  rate  curves  shown  In  Figures  3.12  and  3.13. 

As  for  the  critical  applied  tension  (ax)cr  onset  of  the  mid-plane 
(mode-I)  delamlnatlon,  assume  that  the  effective  Interface  flaw  has  an  Initial 
slae  Uq  -  4t  or  larger.  Then,  the  maximum  values  of  Gg  and  Op  In  Figures  3.12 
and  3.13  are  taken: 


("e)max  =  29.6  x  10*^  J/m^ 
(CT)max  =  J/nVc2 


^  3.5> 


Fro(n  Table  TI-1,  use  aT  =  125°C.,  t  =  0.066  mm  and  Oj-,-;  =  227.5  J/n^. 
Substituting  these  values  Into  Equations  (3.3)  and  (3.4),  It  follows  that  the 
critical  lamlnati  tensile  strain  (-x^cr  onset  of  delamlnatlon  shoull  be 

equal  to  or  greater  than, 

(ex)or  ^  ^-12  x  lO""^  (3.6) 

This  prediction  agrees  with  the  experiments!  Finding  reported  In  Rerereii- 
ces  [32]  and  [41].  Four  replicate  specimens  were  testerl  and  lilamlnatlon 
onset  occiirred  at  laminate  strain  e,^  ranging  Froio  'i.S  to  6.1  x  ]1“^.  The 
average  was  -  5.95  x  10“^. 

In  predicting  (?xX:r  (3.6),  Op,-;  was  nssume.il  to  hao  e  the  coiistMt  yal  ie 
227.5  J/n2.  This  valne  represents  the  upper  l)oun  1  valie  for  ni’-l  in/ii 
Interface  delamlnatlon,  see  Table  ll-I.  For  the  nr>seot  proM'-i,  the  i  ‘  t  i  i  ’ 
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Gic  should  be  somewhat  smaller.  Then,  the  prediction  and  the  experiment  would 
agree  even  closer. 

Quasl-Isotroplc  Laminates  Coupon  Under  Compression.  To  further  Illus¬ 
trate  the  utility  of  the  energy  release  rate  method,  consider  the  following 
laminate  coupons  of  different  lamination  stacking  geometries;  [902/02/^A52lg, 
[O2/9O2/M52IS  and  [ (0/90)2/(+45)2]3*  Note  that  these  laminates  are  macro- 
scoplcally  Identical  In  the  elastic  responses  to  In-plane  laminate  loading 
except  for  the  free  edge  effect.  Specifically,  the  stacking  sequence  of  each 
laminate  causes  tensile  Interlaminar  normal  stress  oz,  shearing  stresses  txz 
and  Tyz  near  the  free  edge  under  uniaxial  compressive  load.  Of  course,  the 
magnitudes  of  the  free  edge  stresses  will  differ  with  the  different  .stacking 
sequences. 

In  Reference  [42],  laminates  of  these  stacking  sequences  were  made  using 
AS-3501-06  graphite-epoxy.  Basic  ply  properties  are  given  In  Table  IIl-l. 
The  laminates  were  tasted  in  room  temperature  under  nonotonlc  uniaxial  com¬ 
pression.  In  all  test  cases,  free  edge  de lamination  In  the  laminate  mid-plane 
was  recorded,  see  Figure  3.4. 

Figure  3.15  shows  the  experimentally  measured  free  edge  delamlnatlon 
growth  plotted  against  the  applied  laminate  compressive  stress  for  each  of 
of  the  three  types  of  laminate  coupons.  These  delamlnatlon  growth  curves 
Indicate  that  onset  of  delamlnatlon  la  followed  by  relatively  rapid  growth, 
and  the  growth  may  have  precipitated  localised  ply  buckling  failure. 

Note  the  onset  load  for  delamlnatlon.  As  can  be  seen  from  Figure  3.15, 
this  varies  with  the  laminate  stacking  sequence.  The  Final  failure  load  Is 
also  affected  by  the  stacking  sequence  factor. 

In  order  to  explain  these  observed  differences,  an  edge  stress  field 
analysis  was  performed  using  the  finite  element  routine.  ^^igures  3.16  -  3.13 
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show  the  Isometric  view  of  distributions  on  the  coupon  cross-section  (first 
quadrant  section)  for,  respectively,  (902/02/±452l s,  (02/902/t452]3  and 
[(0/90)2/(4^5)213.  Similarly,  Figures  3.19  -  3.21  show  the  shearing  stress 
Tx2  fop  same  three  coupon  types.  The  other  shearing  stress,  Is  not 

shown  due  to  Its  relatively  small  magnitude. 

From  these  edge  stress  displays.  It  appears  that  both  the  laminate  mid- 
plane  and  the  A5/-45  Interface  are  possible  delamlnatlon  sites.  It  Is  not 
Immediately  clear,  however,  as  to  why  the  onset  loads  for  delamlnatlon  should 
be  different,  and  why  the  delamlnatlon  sites  shoull  all  at  the  mid-plane. 

To  answer  these  questions,  an  energy  release  rit*  analysis  was  perfi>rmed. 
The  energy  release  rate  coefficient  Function  Go(a)  computed  by  assuming  an 
Interply  effective  flaw  of  slie  a^  la  the  laminate  mid-plane,  and  In  the  45/- 
45  Interface  are  shown  respectively  la  Figures  3.22  and  3,23.  The  Cq—curves 
due  to  residual  thermal  loading  were  neglected  because  of  their  small  values. 

The  mid-plane  delamlnatlon  again  contains  purely  mode-T  action,  while  In 
the  45/-45  Interface  It  contains  both  opening  (mode-l)  and  anti-plane  shear 
(mode-TII)  actions.  The  mixed-mode  ratio  Hm/Gq  varies  from  O.l  to  0.3 
amongst  the  three  types  of  laminates,  see  Reference  [42].  The  Cp-curves  In 
Figure  3.23,  however,  are  for  the  total  energy  release  rate. 

From  the  computed  Cg-ciirves,  It  can  be  Inferred  that  mode-I  (mil-plane) 
delamlantlon  should  occur  In  all  the  three  types  of  laminates.  This  Is 
because  the  computed  Gg-curves  for  mid-plane  delamlnatlon  are  all  higher  than 
the  corresponding  Gg-curves  for  delamlnatlon  In  the  45/-45  Interface  plane. 
In  addition,  for  the  material  system  used,  Gir,  Is  lower  than  mixed-mode 
total,  see  Table  IIl-I. 

The  actual  onset  loads  are  predicted  using  the  maximum  values  of  the 
corresponding  Gg-curves  for  mid-plane  delamlnatlon.  \s  was  done  In  the 
previous  example.  Equations  (3.3)  and  (3.4)  are  used  to  predict  the  minimum 
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possible  onset  strain  e^  for  each  case  considered.  The  results  are  as 
follows: 


Laminate  Type 
[902/02/±452l3 

(02/902/4452!  3 
[(0/90)2/(':45)2ls 


Minimum  Onset 

-303  Mpa  (-44  k,sl) 
-352  Mpa  (-51  ksl) 
-407  Mpa  (-59  ksl) 


The  corresponding  experimental  loads  are  determined  from  Figure  3.15. 
Because  of  the  data  scatter,  load-ranges  are  obtained  yielding,  respectively, 
290  -  310  Mpa;  331  -  358  Mpa;  and  400  -  427  Mpa.  The  predicted  results  are 
within  the  experimental  scatter  ranges.  In  particular,  the  difference  In  the 
onset  loads  caused  by  the  difference  In  stacking  sequences  Is  also  correctly 
predicted. 

In  the  examples  discussed  In  this  section,  the  question  regarding  the 
effective  flaw  size  a^  Is  treated  Indirectly.  The  assumption  used  Is  simply 
that  along  the  free  edge  of  the  laminate,  the  Interface  flaw  has  a  probable 
size  equal  or  greater  than  a,,  (sn.  Is  of  the  order  of  4t  or  the  layer  thickness 
In  which  the  assumed  delamlnatlon  Is  contained).  In  this  manner,  the  predic¬ 
tion  Is  made  using  the  maximum  available  energy  release  rate,  yielding  a  lower 
bound  for  the  delamlnatlon  onset  Toad. 

3 . 4  Two-Olmens Iona 1  Telamlnatlon  Grow th 

In  the  preceding  sections,  the  energy  release  rate  method  has  been 
applied  to  describe  a  number  of  free  edge  delamlnatlon  problems.  The  lami¬ 
nates  studied  were  so  designed  and  loaded  as  tj  Induce  delamlnatlon  without 
significant  Interactions  with  other  cracking  modes.  furthermore,  the  delaml¬ 
natlon  was  either  actually  confined  or  theoretically  assumeil  to  grow  unl'^ornly 


along  the  free  edges  of  the  test  coupon 


80,  It  could  be  treated  as  a  one- 
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dimensional,  self-similar  crack  growth  problem.  Under  such  circumstances,  the 
kinematics  of  the  crack  propagation  Is  simplified  and  crack  growth  Is  repre¬ 
sented  In  a  two-dimensional  plane.  The  crack  size  then  has  a  linear  magnitude 
"a"  and  the  propagation  la  simply  along  a  line.  This  simplification  makes 
possible  a  two-dimensional  stress  field  solution  in  which  the  crack-tip  energy 
release  rate  function  G(^,  a)  Is  computed.  Within  the  general  confines  of  the 
classical  fracture  mechanics,  the  problem  Is  treated,  at  least  conceptually, 
without  recourse  to  empirical  formulisms. 

Generally  speaking,  one-dimensional  delamlnatlon  growth  In  Laminates  Is  a 
very  special  case.  It  can  probably  happen  only  under  Ideally  controlled  condi¬ 
tions.  Most  of  the  delamlnatlon  problems  encountered  In  practice  are  caused 
by  laminate  defects  much  inore  severe  than  free  edges.  The  growth  of  delamlna¬ 
tlon  near  a  local  defect  Is  almost  always  mixed-mode  and  multl-dlrectlonal  In 
nature. 

Vor  example,  consider  the  laminate  with  a  through-hole  shown  In  Plgure 
3.5.  Upon  loading,  localized  delamlnatlon  and  other  modes  of  matrix  cracking 
are  Induced  near  the  curvilinear  edge  of  the  hole.  In  this  case,  the  result¬ 
ing  delamlnatlon  crack  front  Is  actually  a  plane  contour.  Because  the  associ¬ 
ated  growth  Is  two-dimensional,  the  progression  of  the  crack  must  he  represen¬ 
ted  both  In  terms  of  the  Instantaneous  delamlnatlon  area  and  the  delaminated 
shape. 

Clearly,  this  problem  not  only  requires  a  three  dimensional  stress 
analysis  but  also  a  two-dimensional  fracture  growth  treatnent.  Among  the  mist 
difficult  tasks  analytically  are  (I)  how  to  calculate  the  crack  front  euergv 
release  rate  G(o,  a).  In  which  a  Is  some  line  contsur;  and  (2)  now  ti  set  in  a 
fracture  growth  criterion  appropriate  For  the  las t an t  ineous  mixed-mode  ml 
multl-dlrectlonal  crack  growth  hehavlor. 
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If  the  line  contour  of  a  local  de lamination  Is  prescribed  and  Its  Instan¬ 
taneous  growth  shape  assumed,  the  crack  front  energy  release  rate  n(o,  a)  can 
be  calculated  by  a  three-dimensional  finite  element  routine  [28],  Of  course, 
the  calculated  G(o,  a)  remains  as  an  effective,  or  averaged  quantity  In  the 
same  context  as  discussed  In  Chapter  1,  Section  1.5.  As  for  a  criterion 
suitable  for  contoured  crack  propagation  In  a  given  plane,  little  Information 
Is  available  at  the  present  time.  So,  rather  than  attempt  a  rigorous  treat¬ 
ment  encompassing  all  these  fundamental  questions,  consider  the  following 
example  problem  In  light  of  the  analytical  developments  thus  far. 

An  Example  for  Contoured  Delamlnatlon.  Earlier  In  Section  3.3,  free  edge 
delamlnatlons  In  three  similar  quasl-lsotrlplc  laminates  Induced  by  axial 
compression  were  examined.  In  a  related  but  separate  experiment  (421,  one  of 
the  laminate  coupons,  (02/902/14521 st  fabricated  with  a  pre-lmplanted 

delamlnatlon  In  Its  mid-plane,  see  Inset  In  Figure  3,24.  The  experimental 
details  regarding  the  fabrication  and  test  procedures  were  reported  In  Refer¬ 
ence  (42),  Figure  3.24  shows  the  test  results  In  terms  of  the  delamlnatlon 
sl'^e  versus  the  applied  laminate  compression  for  specimens  with  and  without 
the  Implanted  delamlnatlon.  These  laminates  display  strikingly  different 
growth  behaviors.  Specifically,  the  laminates  with  the  Implant  underwent 
contoured  delamlnatlon  emanating  from  the  corner  of  the  Implant  edge,  while 
the  laminate  without  the  Implant  experienced  more  or  less  uniform,  one- 
dlmenslonal  free  edge  delamlnatlon  as  discussed  before.  The  contnnrid  growth 
In  the  laminate  with  Implant  was  comparstlvely  slower  and  more  stshle;  hut  Its 
early  Initiation  and  subsequent  growth  actually  led  to  premature  laminate 
final  failure. 

In  what  follows,  a  prediction  will  be  attempted  for  the  onset  if  the 
contoured  delamlnatlon  growth  with  the  all  <if  the  3-h  finite  element  slmnli- 
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tion  routine  developed  In  [28]. 

But  before  proceeding  with  the  growth  simulation,  It  Is  useful  to  examine 
the  stress  field  near  the  Implant  at  the  coupon's  free  edge.  For  In  this 
region  a  severe  stress  concentration  Is  present.  The  particular  point  of 
Interest  Is  the  Interface  corner  where  the  Implanted  delamlnatlon  intersects 
the  free  edge  of  the  test  coupon.  This  point  Is  labeled  "Interface  corner"  In 
Figure  3.25. 

Recall  that  laminates  without  the  Implant  developed  a  tensile  Interlami¬ 
nar  oz  stress  near  the  laminate  free  edge.  The  largest  value  of  oj  near  the 
free  edge  (acting  on  the  mid-plane)  was  about  14  kpa  (2  psl)  per  far-fleld 
strain  e^  =■  -10“^;  see  Figure  3.16.  In  contrast,  for  laminates  with  the 
Implant,  the  distribution  of  oz  on  the  laminate  mtd-plane  Is  shown  In  Figure 
3.25,  where  the  magnitude  of  oz  near  the  "Interface  corner"  Is  In  the  order  of 
60  kpa  (8  psl)  .  per  =•  In  fact,  oz  the  point  Is  theoretically 

unbounded. 

Since  the  Implant  Is  situated  la  the  lamlnte  mid-plane,  the  Initial 
delamlnatlon  growth  Is  primarily  of  mode-I. 

To  predict  the  associated  delamlnatlon  onset.  It  Is  necessary  to  calcu¬ 
late  the  crack-front  strain  energy  release  rate.  Assume  that  the  expected 
Initial  delamlnatlon  will  start  from  the  Interface  corner  point  and  grow  with 
a  certain  known  contoured  front.  This  growth  contour  Is  represented  by  a 
prescribed  sequence  of  nodal  releases  In  the  finite  element  mesh.  Specific¬ 
ally,  Figure  3.26(a)  shows  the  Intended  nodal  release  sequence  that  'dmlcks 
the  delamlnatlon  growth  In  a  discrete  fashion.  The  actual  number  labeled  for 
each  node  represents  the  order  of  its  release.  In  this  manner,  the  energv 
release  rats  coefficient  rig  can  be  computed  as  a  function  of  the  dclaminatad 
area  associated  with  the  prescribed  growth  contaur.  This  Is  shown  la  Figure 
3.26(b). 
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It  Is  seen  that  the  delamlnatton  growth  Is  tnost  energetic  along  the 

Implant  edge,  represented  by  the  nodes  numbered  1,  2,  3,  *>,  6,  8,  etc.  As  the 
delaminated  area  becomes  larger,  the  available  energy  release  rate  becomes 
smaller.  This  Implies  that  the  growth  process  Is  essentially  stable,  as 
observed  In  the  experiment. 

To  obtain  a  reasonable  prediction,  the  maximum  value  of  Cg  In  Figure 
3.26(b)  Is  used  to  calculate  the  onset  load  for  delamlnatton.  In  this  case, 
thermal  effect  on  delamlnatton  Is  Ignored  as  Is  relatively  small.  Then, 
the  maximum  strain  energy  release  rate  associated  with  the  orescrlbed  initial 
delamlnatton  contour  can  be  expressed  as: 

Gmax  “  (Gg),nax  (2t)  (3.7) 

where  the  value  of  (Ge),nax  corresponds  to  the  release  of  node  6  la  Figure 

3.26(b). 

Using  Og  =  175  J/;n^  (1.0  In-lb/ln^)  for  mode-T  delamlnatlon  along  the 

fiber  direction  and  2t  =  0.13  mm  (0.005  In.)  from  Table  Tll-l,  the  calculated 
minimum  onset  load  In  terms  of  the  applied  far-fleld  strain  (e^)  Is  found  to 
be 

(ex)mln  =  2.80  x  10"^  (1.8) 

For  the  laminate  considered,  the  axial  stiffness  Is  =  50.1  Opa  (7.27 
msl),  so  that  the  predicted  (minimum)  onset  far-fleld  laminate  stress  Is 

about 

(cx^mln  "  (21.0  ksO  (3.0) 

This  predicted  value  Is  almost  the  same  as  the  experimental  value  oxtrec- 
ted  from  Figure  3.25. 
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Up  to  this  point,  reasonable  analytical  rigor  could  be  nialntalned.  Hut, 
what  to  do  when  the  delamlnatlon  growth  becomes  larger  Is  not  at  all  clear. 
Hot  only  Is  the  exact  crack  front  unknown  but  the  growth  also  becomes  mixed¬ 
mode. 


I.**' 

‘'i: 


3.5  Summary 

In  this  chapter,  the  basic  mechanics  of  Interply  cracking  has  been 


discussed  from  both  an  experimental  and  an  analytical  view  point.  Several 
examples  belonging  to  the  class  of  free  edge  Induced  de laminations  were 
studied  by  means  of  the  energy  release  rate  method  In  conjunction  with  a  three 
dimensional  ply  stress  analysis.  On  the  whole,  the  energy  method  Is  a  useful 
technique  to  describe  some  of  the  essential  characteristics  found  experiment¬ 
ally  In  free  edge  delamina tlons. 

The  basic  mechanisms  of  Interply  cracking  seem  to  be  more  Intricate  than 


Intraply  cracking.  Since  an  Interply  crack  Is  propagating  parallel  (albeit 
macroscoplcally)  to  the  laminate  plane,  the  applied  In-plane  load  required  to 
propagate  an  Interply  crack  Is  usually  higher  than  that  required  to  propagate 
an  Intraply  crack.  For  this  reason,  Interply  cracking  almost  always  Interacts 
with  other  matrix  cracking  events  In  practical  cases. 

Even  In  the  special  case  of  free-edge  delamlnatlon,  such  as  those 
examples  studied  In  this  chapter,  one  cannot  always  create  a  pure  delamlnatlon 
growth  without  concurrent  Interactions  with  other  cracking  modes.  To  model 
such  Interactions  would  necessarily  add  uncertainties  both  In  physical  reason¬ 
ing  and  In  analytical  modeling. 

An  equally  troublesome  aspect  In  the  delamlnaton  problem  Is  the  lick  of  a 
suitable  mixed-mode  growt'i  criterion.  Because  of  the  anisotropic  natiri  of 
the  composite  piles,  an  Interface  crack  woul1  always  propagate  In  ml<ed-mode. 


And,  for  the  same  reason,  the  material  resistance  against  crack  propagation 
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also  depends  on  the  ply  geometry  (e.g.  fiber  orientations)  as  well  as  the 
particular  mlcrostructures  that  exist  on  the  Interface. 

<111  these  difficulties  make  It  seemingly  Impossible  to  treat  the  problem 
on  a  more  rigorous  basis.  Nonetheless,  further  analytical  developments  remain 
promising.  \s  will  be  demonstrated  In  the  next  chapter,  the  energy- 
release-rata  method,  as  outlined  here,  continues  to  provide  a  reasonable 
analytical  description  for  matrix  crack  Interactions  that  emerge  at  higher 
applied  load  levels. 
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TABLE  TTI-1 

Material  Constanta  Cor  AS-3501-06  Untdtrecttonal  Ply* 


Property 

SI  Unit 

English  Unit 

^LL 

140.0  Gpa 

2O.4x106  psl 

^TT* 

ll.O  Gpa 

l.hxlO^  psl 

^Ti  VT,z 

0.29 

0.29 

'Tz 

0.3 

0.3 

GlT«  “"’Lz 

5.5  Gpa 

O.BxlO^  psl 

GTz 

5 . 5  Gpa 

0.8x10^  psl 

QT»  az 

28.Bx10"6/°C 

16.l0-^/°F 

OL 

o.ie.io-^/^c 

0.2xl0”'’/°F 

OL,  long,  tensile  strength 

1325  mpa 

190  ksl 

OT,  trans.  tensile  strength 

50  mpa 

3.5  ksl 

2t,  nominal  ply  thickness 

0.132  mm 

0.0052" 

at,  temperiture  drop 

155'’C 

3l0'’F 

Gjc  (G/0  Interface)** 

175  J/.n^ 

l.O  In-lh/ln- 

'jc  (90/90  Interface)** 

223  I/n2 

1.3  in-lb/ln^ 

♦Data  taken  from  fere ace  [42], 
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X-Radiograph  (left)  and  Micrograph  Showing  a  Cross-Sectional 
Delamination  in  a  [+25/90/+25]  Tensile  Coupon. 
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I  Figure  3.5  X-Radiograph  Showing  Matrix  Cracking  Around 
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Figure  3.11  Through-Thickness  Distribution  of  Near  the  Free  Edge  of  a  [+25/90] 
Tensile  Coupon. 
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Computed  Energy  Release  Rate  Coefficient  Cg(a)  for  the  [+25/90]g  Tensile 


Figure  3.13  Computed  Energy  Release  Rate  Coefficient  Cj(a)  for  the  [+25/90]g  Tensile 
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Figure  3.15  Experimental  Free  Edge  Delamination  Growth  Plotted  Against  the  Applied 
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Figure  3.22  Mode  I  Energy  Release  Rate  Coefficient  C^Ca)  for  all  Three  Compression  Coupons. 
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Figure  3.26  The  Proposed  Initial  Delamination  Nodal  Release 
Sequence  (above)  and  the  Computed  Energy  Release 
Rate  Coefficient  Cg(A). 
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CHAPTER  IV.  ANALYSIS  OF  LATE-STAGE  MATRIX  CRACKING 

4.1  Matrix  Cracks  at  High  Load  Levels 

From  the  analysis  In  the  previous  chapters,  It  may  be  concluded  that 
matrix  damage  begins  generally  with  the  formation  of  Intraply  cracks  In  the 
secondary  or  non-loading  plies.  In  one  example  discussed  in  Chapter  TI,  the 
cross-ply  l02/904]s  laminate  suffered  sporadic  transverse  cracks  In  the  h0°- 
plles  at  laminate  strains  as  low  as  0.2%  to  0.3%.  Since  many  laminates 
designed  for  practical  usage  can  reach  ultimate  tensile  strains  In  excess  of 
1%,  early  stage  intraply  cracking  In  the  secondary  plies  Is  often  tolerated  as 
benign  service-related  defects. 

In  some  Instances,  interlaminar  matrix  cracking  In  the  form  of  Jelamlna- 
tlon  can  occur  before  Intraply  cracking.  Such  occurrences  are  mostly  confined 
to  localized  areas,  such  as  near  the  laminate  free  edges,  around  man-made 
notches  or  at  the  leading  edges  of  Internal  damages,  where  high  concentration 
of  Interlaminar  stresses  Is  present. 

In  any  event,  both  types  of  matrix  cracking  are  undesirable  from  the 
laminate  reliability  and  durability  viewpoint.  This  Is  because  matrix  cracks 
In  secondary  piles  can  precipitate  matrix  cracks  Into  the  load-carrying  olles, 
under  either  monotonlcally  Increasing  or  cyclically  applied  loads. 

In  the  cross-ply  laminates  discussed  l.a  Chapter  IT,  for  instance,  trans¬ 
verse  cracks  In  the  h0°-plles  did  not  seem  to  affect  the  structural  perform- 
ance  of  the  laminate  until  about  h0%  the  ultimate  loail  [38].  Toward  the  late- 
stage  of  loading,  however,  In-ply  matrix  splits  In  the  n°-plles  were  found 
throughout  the  Interior  of  the  laminates.  The  splits  are  orientated  parillel 
to  the  applied  tension  and  cross-over  the  'i0°-ply  transverse  cracks  that 
formed  early.  4s  soon  as  a  n°-ply  split  formed  and  crossed  over  a  h'^-nlv 
transverse  crack,  localized  delamlnatlon  in  the  0/')0  Interface  was  found  near 
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the  crossing.  Subsequent  rapid  propagation  of  many  of  the  cross-cracks 
Induced  delamlna tlons  eventuated  a  massive  0/90  Interface  separation.  This 
event  was  followed  almost  Immediately  ay  catastrophic  failure  of  the  0®-plles. 
The  last  failure  event  Involved  fiber  breaking  In  the  split  0®-plles.  Figure 
4.1  shows  a  sequence  of  x-ray  photographs  that  Illustrate  the  matrix  cracking 
proliferation  process  of  a  [02/90313  laminate  during  Its  late-stage  loading. 
These  photographs  clearly  demonstrate  that  matrix  cracks  In  non-loading  piles 
can  become  detrimental  to  the  load-carrying  piles  at  high  load  levels. 

Laminates  of  more  practical  lay-ups  suffer  similar  late-stage  matrix 
cracking,  though  the  sequence  of  events  may  not  be  the  same  [47].  Experiments 
using  small  coupons  confirmed  that  late-stage  matrix  cracking  patterns  are 
usually  generic  to  the  lamination  geometry  (stacking  sequence,  ply  thickness, 
etc.)  and  the  nature  of  loading  [47,  48].  The  same  factors  that  are  Important 
In  the  early-stage  matrix  cracking  continue  to  play  a  dominate  role  during  the 
late-stage  matrix  cracking. 

The  objective  of  this  chapter  Is  to  Investigate  some  of  the  basic  mechan¬ 
isms  associated  with  late-stage  matrix  cracking.  The  basic  macromechanlcs/ply 
elasticity  approach  discussed  In  the  previous  chapters  is  followed.  Since  the 
cracks  concerned  Involve  three  dimensional  and  localised  stress  fields,  the 
finite  element  routine  developed  In  Reference  [28]  Is  employed  In  the 
analysis. 

In  order  to  achieve  a  cohesive  discussion,  emphaslTilng  mechanisms  rather 
than  the  analysis  method,  a  specific  family  of  laminates  will  be  analyzed. 
Namely,  the  T300/934  tensile  coupons  4  are  Investigated 

again  for  late-stage  matrix  cracking.  A  complete  documentation  of  damage 
accumulation  up  to  final  laminate  failure  Is  available  In  Reference  f38l;  so  a 
comparison  between  the  analysis  and  experiment  can  he  made. 
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4.2  Stress  Field  Near  a  Transverse  Crack  Terminus 

Consider  again  the  [02/90^13  laminate  coupon  series  under  monotonlc  axial 
tension.  Transverse  cracks  In  the  90®-plle8  are  formed  during  the  early-stage 
loading.  The  effect  of  these  cracks  on  damage  development  during  the  late- 
stage  loading  can  be  qualitatively  studied  by  examining  the  behavior  of  the 
stress  fields  near  the  transverse  crack  termini. 

Figure  4.2  shows  a  tensile  coupon  which  already  suffered  90‘’-ply  trans¬ 
verse  cracks.  For  the  purpose  of  presentation  only,  consider  the  section  of 
the  laminate  that  contains  only  one  transverse  crack.  Due  to  lamination 
symmetry,  only  one-eighth  of  this  section  needs  to  be  modeled,  see  Figure  4.2. 
Details  of  the  finite  element  sizing,  mesh  generation  and  other  computational 
aspects  are  found  In  Reference  (28]  and  will  not  be  discussed  here. 

The  computed  stress  field  near  the  transverse  crack  terminus  Is  complex, 
but  the  stresses  can  still  be  separated  In  two  groups.  The  first  group 
consists  of  the  Interlaminar  stresses  02,  t^z  tyz  acting  on  the  0/90 
Interface.  These  stresses  are  responsible  for  delamlnatlon.  The  second  group 
consists  of  the  In-ply  stresses  o^t  oy  and  i^y  acting  In  the  0“-plles.  These 
stresses  may  cause  matrix  damage  In  the  O^-ply,  such  as  length-wise  splitting, 
as  well  as  fiber  breaking. 

Interlaminar  Stresses.  Figures  4.3  to  4.5  show,  respectively,  the  Inter¬ 
laminar  Oz,  T,cz  fyz  distributions  on  the  0/90  Interface.  In  each  display, 
the  x-y  plane  coincides  with  the  0/90  interface;  the  x-axls  coincides  with  the 
free-adge  of  the  coupon,  while  the  loci  of  the  transverse  crack  terminus  Fall 
on  the  y-axls.  One  of  the  corner  points  labeled  In  Figure  4.2  Is  located  at 
the  x-y  origin. 

In  each  display,  stresses  for  the  l32'^^^2l  s  ■''’^d  [02/90/,]  3  coupons  are 
shown.  The  magnitude  of  the  stresses  Is  due  to  an  applied  laminate  strata 

®x  ^  1  pa.  It  is  noted  that  the  Interlaminar  stress  fields  for  a  =  2  and 


a  ••  4  show  similar  distributional  features,  but  with  noticeable  difference  In 
magnitudes  and  In  areas  of  concentration.  In  particular,  the  Interlaminar 
normal  stress  oz  Is  tensile  In  both  cases,  concentrated  along  the  x  and  the  y 
axes.  This  peculiar  distribution  Is  effected  by  both  the  free  edge  and  the 
transverse  crack  In  the  QO^-layer. 

The  tensile  oz  and  the  shearing  stresses  txz  and  lyz  could  cause  delaml- 
natlon  In  the  0/90  Interface,  emanating  from  the  corner  point  and  growing  Into 
a  delta-shaped  contour  as  Illustrated  In  Figure  4.6.  \n  analysis  of  that 
problem  will  be  presented  In  the  next  section. 

In-Plane  Stresses  In  O^-Plles.  As  for  the  tn-ply  stress  group  -  ox*  Oy 
and  Txy,  especially  those  existing  In  the  0°-plies,  their  intensities  are  also 
affected  by  the  presence  of  the  free  edge  and  the  transverse  crack. 

First,  consider  the  distributional  characteristics  of  these  stresses  In  a 
plane  parallel  to  the  0/90  Interface,  or  the  x-y  plane.  Figures  4.7  to  4.9 
show,  respectively,  Oxt  ay  Txy  distributions  In  the  0®-ply  near  the  0/90 
Interface.  A  severe  stress  concentration  occurs  near  the  transverse  crack 
terminus  and  a  less  severe  concentration  occurs  along  the  free  edge  la  all 
cases.  The  90°-layer  thickness,  or  the  value  of  n,  is  seen  to  Influence  both 
the  stress  concentration  Intensity  and  the  size  of  the  stress  concentration 
zone. 

Of  the  three  In-plane  stresses,  the  ■nagnltude  of  ixy  Figure  4.9  is 
relatively  small,  though  Its  illstrlbutlonnl  behavior  Is  quite  dramatic  near 
the  corner  point.  However,  ox  and  especially  oy,  magnified  significantly  by 
the  presence  of  the  transverse  crack,  will  cause  damage  In  the  0°-ply.  In 
particular,  both  stresses  are  In  tension  throughout  the  flell  and  each  can 
causa  a  different  mode  of  failure  of  the  0°-ply. 

The  thernal  resHual  stresses  due  to  thermal  mismatching  of  the  0°  and 
90“  piles  are  also  significant.  For  reasons  of  brevity,  they  are  not  shown. 
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but  their  distributional  features  are  similar  to  those  displayed  In  Figures 
4.3  to  4.5,  and  4.7  to  4.9.  When  coupled  with  the  applied  load,  the  role  of 
the  thermal  stresses  In  matrix  damage  development  cannot  be  Ignored. 

The  distributional  behavior  of  ox  ^nd  oy  In  the  0®-layer  near  the  trans¬ 
verse  crack  terminus  cannot  be  fully  revealed  by  displays  such  as  Figures  4.7 
and  4.8.  Both  stresses  also  vary  dramatically  throughout  the  thickness  of  the 
O^-layer. 

Figure  4.10  shows  the  Ox  distribution  In  the  0°-layer  thickness  above  the 
transverse  crack  terminus  (see  figure  Inset).  Note  that  the  part  due  to  the 
applied  tension,  ex  =  1  |ie,  Is  tensile  throughout  the  thickness  and  Is  highly 
concentrated  approaching  the  0/90  Interface.  The  part  due  to  thermal  cooling, 
AT  =  “l“f,  Is  compressive,  except  near  the  transverse  crack  root  region  where 
It  becomes  tensile  with  significant  tnagnltude.  The  combined  ox  near  the 
transverse  crack  terminus  region  displays  varied  degree  of  concentration, 
depending  on  the  thickness  of  the  00®-layer  In  which  the  transverse  crack  is 
situated. 

It  Is  apparent  that  concentrations  of  Ox  nnnr  the  transverse  crack  sites 
may  have  a  profound  effect  on  the  breaking  strength  of  the  0°-ply  fibers. 
Several  recent  experiments  [47,  4Sl  found  broken  fibers  in  the  0“-plies  at 
locations  where  adjacent  OO^-ply  transverse  cracks  had  occurred. 

Figure  4.11  displays  a  similar  distribution  for  Oy,  which  is  an  in-plane 
stress  normal  to  the  fibers  of  the  0°-ply.  It  is  seen  that  the  part  of  oy  due 
to  Ox  =  I  pe  and  the  part  due  to  AT  =  -1°F  are  both  tensile  throughout  the 
thickness  of  the  0°-layer.  Both  are  magnified  la  intensity  near  the  9n°-ply 
transverse  crack  terminus.  The  near-fleH  oy  distribution  is  seen  to  '’ary 
again  with  the  OQ^-ply  thickness  or  the  value  of  n.  The  high  tensile  stress 
is  capable  of  causing  matrix  cracks  in  the  n'’-layer  in  the  form  o^  longitndl- 
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nal  splits.  The  splits  occur  in  much  the  same  manner  as  the  90'*-Dly  trans¬ 
verse  cracks  discussed  in  Chapter  II. 


4,3  Damage  Modes  Effected  by  Transverse  Cracks 

From  the  preceding  section,  it  may  be  summarized  that  the  presence  of 
transverse  cracks  in  the  [02/^OyjIs  laminates  can  cause  three  possible  damage 
modes,  all  of  which  are  localized  near  the  transverse  crack  termini.  Of  the 
three  possible  modes,  the  0/90  Interface  delamination  and  the  0°-ply  longitu¬ 
dinal  splitting  are  matrix-dominated  cracks,  while  C’-ply  fiber  breaking  is 
something  different  fundamentally.  The  subject  of  fiber  breaking  requires  a 
different  theoretical  treatment  and  is  outside  the  scope  of  this  discussion. 
Hence,  only  the  localized  delamination  and  0“-ply  splitting  are  analyzed  in 
this  section. 

Onset  of  Delamination  at  Transverse  Crack  Terminus.  The  physical  reason 
for  delamination  at  a  transverse  crack  terminus  Is  quite  clear.  The  malor 
stresses  causing  delamination  in  the  0/90  interface  are  the  concentrated 
tensile  and  shearing  xxz*  other  shear  stress  ryz  appears  to  be  negli¬ 

gible,  see  Figures  4.3  to  4.5.  These  stresses  attain  their  highest  magnitude 
at  the  corner  point  where  interactions  take  place  between  the  transverse  crack 
and  the  free  edge,  see  Figure  4.6.  Thus,  if  a  delamination  is  to  occur,  it 
would  probably  start  from  this  point  and  propagate  into  a  plane  contour  as 
shown  in  the  figure. 

In  order  to  make  a  prediction  for  the  onsat  and  the  growth  behavior  of 
the  contoured  delamination,  a  small  Interface  effective  flaw  will  again  be 
assumed,  located  near  the  corner  point.  This  assumption  is  in  accordance  with 
the  effective  flaw  concept  discussed  earlier.  Then,  under  the  applied 
tension,  the  assumed  flaw  starts  to  propagate  and  becomes  a  local  delamina¬ 
tion.  Following  the  energy  release  rate  methodology  outlined  before,  the 

no 
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growth  of  the  crack  front  can  be  mimicked  by  a  sequence  of  nodal-point 
releases  (Figure  4.6)  In  the  finite  element  net-work,  shown  In  Figure  4.2.  4s 
each  finite  element  node  Is  released,  the  corresponding  delamlnatlon  area  as 
well  as  the  strain  energy  lost  can  be  approximately  computed  (for  details,  see 
References  [28]).  In  this  manner,  a  relationship  between  the  delaminated  area 
and  the  strain  energy  release  rate  Is  obtained  for  the  assumed  delamlnatlon 
growth. 

Figure  4.12  shows  the  computed  energy  release  rate  coefficient  Cg  (for 
ex  =*  1)  plotted  against  the  assumed  delamlnatlon  area  contour  for  the 
[O2/9O2IS  laminate.  Note  that  Cg  consists  of  all  three  modal  components  of 
cracking.  In  this  case,  the  mode-I  (z)  component  Is  slightly  larger  than  the 
mode-II  (y)  component,  while  the  mode-IIl  (x)  component  Is  negligible. 

Figure  4.13  shows  the  Ct  coefficient  for  the  laminate  subjected  to 
at  *  -l^F  (or  -0.65®C).  The  general  character  of  and  Its  three  modal 
components  of  cracking  are  similar  to  those  of  Cg.  The  calculated  Cg  and  Cp 
coefficients  are  based  on  the  assumed  nodal  release  sequence  (shown  la  Figure 
4.6),  which  may  not  be  the  best  representation  for  the  actual  delamlnatlon 
growth.  But,  the  general  trend  of  these  coefficient  functions  Indicates  that 
the  Initial  delamlnatlon  growth  Is  basically  stable,  as  Is  the  case  obser¬ 
ved  In  experiment. 

Now,  assume  that  the  starting  size  of  macroscoplcally  observable  delaml¬ 
natlon  Is  about  4t^  (where  2t  Is  one  ply  thickness)  so  that  when  such  a  delam¬ 
lnatlon  Is  Induced,  It  defines  the  onset  of  the  contoured  delamlnatlon.  Then, 
by  taking  the  values  of  Cg  and  Op  from  Figures  2.12  and  2.13  at  4/t^  *  4,  the 
total  energy  release  rate  C(ex,  AT)  can  be  computed  via  Equation  (3.3).  bet, 
as  before,  AT  =  12S®C  and  Cg  =  289  J/n^-  for  mixed-mode  crack  growth,  see  Table 
II-l.  It  follows  from  Equation  (3.4)  that  the  predicted  threshold  tensile 
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strain  foe  the  onset  of  the  contoured  delamlnatlon  for  n  >■  2  Is: 

(ex)de  “  11*0  *  10”^ 

Similarly,  the  computed  (Jg  and  Cj»  coefficient  functions  for  [02/90413  are 
shovm  In  Figures  4.14  and  4.15,  respectively.  It  Is  seen  that  the  values  of 
Cg  and  Of  are  significantly  Increased  as  a  result  of  doubling  the  number  of 
the  90“-plle8.  Consequently,  the  predicted  threshold  strain  for  onset  of  the 
assumed  contoured  delamlnatlon  In  this  laminate  Is: 

(ex)de  =“ 

Since  the  breaking  strain  of  the  0®-ply  Is  normally  at  about  10  x  10“^, 
It  may  be  Inferred  that  the  assumed  contour  delamlnatlon  would  not  occur  In 
the  t02/902]s  laminate,  but  vjould  occur  In  the  [02/90413  laminate. 

Indeed,  In  the  experiment  reported  In  Reference  [381,  no  such  delamlna¬ 
tlon  was  detected  until  after  laminate  failure  In  [02/90lg  and  [02/90218* 
[O2/9O4I3,  however,  delta-shaped  delamlnatlon  were  found  at  laminate  strain 
about  7.5  X  10”^.  This  compares  closely  with  the  predicted  value  In  (4.2). 

Onset  of  0°-Ply  Splitting.  As  was  discussed  In  the  previous  section,  the 
basic  mechanisms  In  O^-ply  splitting  are  much  the  same  as  In  90®-ply  trans¬ 
verse  cracking.  However,  the  norma l-to-the-flber  tension,  oy.  In  the  0®-ply 
Is  more  concentrated  near  the  transverse  crack  terminus.  In  addition,  along 
the  length  of  the  laminate  such  concentrations  occur  periodically,  coinciding 
with  the  regularity  of  the  90®-ply  transverse  cracks.  This  makes  It  difficult 
to  predict  how  a  0®-ply  split  Initiates. 

Let  It  be  assumed  that  an  effective  In-ply  flaw  In  the  0“-ply  exists  and 
lies  fiber-wise  In  the  ply.  Then,  the  growth  conditions  for  the  effective 
flaw  can  be  formulated  based  on  the  same  concept  of  energy  release  rate.  48 
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In  the  case  of  the  90“-ply  transverse  cracking,  the  size  of  the  assumed  effec¬ 
tive  flaw  Is  not  known.  But,  by  computing  the  maximum  possible  energy  release 
rate  available  along  the  growth  path  of  such  a  flaw,  a  minimum  possible  onset 
load  for  the  flaw  growth  can  be  obtained.  The  actual  computational  procedure 
Is  Identical  to  the  transverse  cracking  cases  discussed  In  Chapter  II;  there¬ 
fore,  the  details  of  the  computation  will  be  omitted  here.  Instead,  only  the 
final  results  will  be  given.  Namely,  longitudinal  splits  In  the  Interior  of 
the  O^-ply  will  probably  occur  over  the  transverse  cracks  when  the  applied 
tensile  strain  e^  exceeds  the  following  critical  values: 


for  n  =  2, 

(ex)3p  =  10.5  X  10-3 

(4.3) 

for  n  =  4, 

(ex)3p  =  7.8  X  10-3 

(4.4) 

According  to  this  prediction,  0®-ply  splitting  cannot  occur  In  the 
[02/902]g  before  laminate  failure;  but  may  occur  In  the  192^^0413  l^^nlnate 
because  the  approximate  threshold  breaking  strain  for  0“-ply  fiber  Is  about 
10  X  10“3. 

Note  that  the  threshold  strain  for  O^-ply  splitting  In  the  [02/904]g 
laminate  Is  about  the  same  as  the  threshold  strain  for  the  0/90  delamlnatlons 
Induced  near  the  transverse  crack/ free-edge  corner  points,  see  Equation  (4.?). 
Or,  these  two  types  of  matrix  cracks  may  occur  at  the  same  load  level.  This 
was  actually  observed  In  the  [02/903lg  laminates  shown  In  Figures  3.6  and  4.1. 
In  Figure  3.6,  local  delamlnatlons  occurred  near  the  Free  edge,  while  la 
Figure  4,1  0°-ply  splits  occurred  la  the  Interior;  both  types  of  matrix  crack¬ 
ing  occurred  at  about  the  same  applied  load  level.  In  the  experiment  reported 
In  [38],  longitudinal  splitting  was  also  found  In  [O2/9O2I3  at  the  limfnate 
strain  of  9.2  x  10"^,  while  la  (02/90413  at  the  laminate  strains  of  less  than 
B  X  10”^.  These  numbers  are  In  general  agreement  with  the  predictions  In 


Equations  (4.3)  and  (4.4). 
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4.4  Delaminatlon  at  Two  Cross  Inc;  Cracks 


When  a  splitting  crack  is  formed  in  the  0“-ply,  it  makes  a  cross  with  the 
transverse  crack  in  the  90“-ply  as  illustrated  by  Figure  4.16,  In  the  close 
vicinity  of  the  crossing  point,  a  severe  stress  concentration  is  developed. 
Both  interlaminar  and  In-plane  stresses  are  present.  The  interlaminar 
stresses  in  particular  will  cause  localized  delamination  such  as  shown  by  the 
shaded  area  in  Figure  4.16.  (See  also  Figure  4.1.) 

Figure  4.17(a)  shows  the  0/90  interface  plane,  with  the  x-axls  parallel 
to  the  0“-ply  split  and  the  y-axls  parallel  to  the  90®-ply  transverse  crack. 
The  axis  origin  is  the  intersection  point  of  the  two  crossing  cracks. 

The  stress  field  near  the  intersection  point  Is  highly  amplified.  In 
particular,  the  interlaminar  normal  stress  0%  Is  tensile  due  both  to  the 
applied  tension  and  thermal  cooling. 

To  determine  the  threshold  condition  for  delaminatlon  caused  primarily  by 
the  interlaminar  stresses,  a  crack  growth  simulation  is  performed  using  the 
same  finite  element  nodal  release  procedure.  The  simulated  delaminatlon 
follows  the  nodal  release  sequence  shown  in  Figure  4.17(b).  The  polnt- 
by-polnt  energy  release  rate  coefficients  Cg  and  Cj  are  shown  in  Figures  4.18 
and  4.19  for  (02/'^^2l3  snd  in  Figures  4.70  and  4.71  for  respective¬ 
ly. 

From  the  energy  release  rate  coefficient  curves,  it  is  apparent  that  the 
behavior  of  the  delaminatlon  la  mixed-mode.  The  total  energy  release  rite 
slightly  decreases  with  the  delaminated  area  indicating  a  stable  growth 
behavior  at  least  initially. 

Comparing  the  four  energy  release  rate  curves  with  those  shown  la  Figures 
4.12  to  4.15,  a  close  similarity  is  seen  between  the  delaminatlon  occurring  it 
the  ciosslng  cracks  and  at  the  transverse  crack/free-edge  corner  points. 
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However,  the  calculated  energy  release  rate  for  delamlnatlon  at  the  cross 
cracks  Is  slightly  higher. 

Using  the  same  procedure  as  before  for  the  threshold  strain  at  onset  of 
delamlnatlon,  the  predicted  laminate  strains  are  as  follows: 


for  n  “  2 

(ex)de  =■  X  10'^ 

(4.5) 

for  n  =  4 

(ex)de  =  7.5  X  10-3 

(4.6) 

This  shows  that  delamlnatlon  at  the  cross  cracks  occurs  at  about  the  same 
load  as  the  0®-ply  splitting.  Experimentally,  however,  the  splitting  occurred 
first,  while  the  delamlnatlon  followed  Immediately  [381.  Since  the  analysis 
Involved  considerable  approximations,  both  In  numerical  computation  and  In 
setting  the  threshold  conditions  for  flaw  growth,  the  close  agreement  achieved 
between  the  prediction  and  experiment  was  rather  surprising. 

From  these  results,  the  energy  release  rate  concept  coupled  with  the 
effective  flaw  postulate  Is  shown  to  be  useful  In  predicting  a  threshold  load 
for  matrix  cracking  that  takes  place  at  high  loads. 

4. 5  Effects  of  Matrix  Cracking  on  laminate  Strength 

As  has  been  shown,  matrix  cracks  and  crack  Interactions  In  laminates  can 
occur  In  multiple  piles  under  high  loads.  These  cracks  may  severely  reduce 
laminate  strength  if  the  laminate's  Internal  microstructure  (e.g.  ply  stacking 
sequence)  allows  them  to  grow,  either  by  self-propaga tlon  or  by  mutual  coale¬ 
scence.  For  example,  matrix  cracks  In  off-axis  piles,  such  as  hO^-ply  trans¬ 
verse  cracks,  are  Initial  laminate  damages.  In  time,  or  at  high  loads,  these 
cracks  cause  secondary  matrix  cracks  In  the  load-carrying  piles  which  trigger 
catastrophic  laminate  failure.  Since  matrix  cracking  Is  Influenced  greatly  by 
the  lamination  geometry,  material  system  and  loading  environment,  even  a 
qualitative  understanding  of  their  Inter-rela tlonshlps  Is  useful  for  optimum 
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laminate  design.  It  Is  shown  here  that  these  factors  can  be  explicitly 
accounted  for  In  the  predictive  models.  In  the  examples  presented,  It  seems 
that  the  energy  method  can  be  an  effective  means  to  analyze  certain  matrix 
cracking  processes  and  their  contribution  to  the  final  failure  of  the  load¬ 
carrying  piles. 

Up  to  now,  most  practical  design  analysts  have  Ignored  the  detailed 


matrix  cracking  patterns  In  the  off-axis  piles.  They  assume  that  all  loads 
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Figure  4.5  Xyz  Distribution  on  0/90  Interface  Near  a  Transverse 
Crack. 


141 


KPo 

PS! 

hi 

“./6 

.83 

1 

.55 

-.0  8 

.P8 

-.04 

0 

-  0 

Figure  4 


ea 


NADC-85118-60 


Figure  4.12  Strain  Energy  Release  Rate  Coefficient,  Cg,  for  Delamination  at  Corner 
Point  [02/902]s- 
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Figure  4.13  Strain  Energy  Release  Rate  Coefficient,  Cj,  for  Delamination  at  Corner 
Point  [02/902]g. 


.15  Strain  Energy  Release  Rate  Coefficient,  Ct,  for  Delamination  at  Corner 
Point  [O2/9O4]  . 


Strain  Energy  Release  Rate  Coefficient,  C^,  for  Delaminat ion  at  Cross-Cracks 


Release  Rate  Coefficient,  Ct,  for  Delamination  at  Cross-Cracks 
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Figure  4.20  Strain  Energy  Release  Rate  Coefficient,  C  ,  for  Delamination  at  Cross-Cracks 
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Strain  Energy  Release  Rate  Coefficient,  C^,  for  Delamlnatlon  at  Cross-Cracks 
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